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This review presents a summary and evaluations of the superconducting properties of the layered ruthen- 
ate Sr2Ru04 as they are known in the autumn of 2011. This paper appends the main progress that has been 
made since the preceding review by Mackenzie and Maeno was published in 2003. Here, special focus is 
placed on the critical evaluation of the spin-triplet, odd-parity pairing scenario applied to Sr2Ru04. After 
an introduction to superconductors with possible odd-parity pairing, accumulated evidence for the pair- 
ing symmetry of Sr2Ru04 is examined. Then, significant recent progress on the theoretical approaches to 
the superconducting pairing by Coulomb repulsion is reviewed. A section is devoted to some experimental 
properties of Sr2Ru04 that seem to defy simple explanations in terms of currently available spin-triplet sce- 
nario. The next section deals with some new developments using eutectic boundaries and micro-crystals, 
which reveals novel superconducting phenomena related to chiral edge states, odd-frequency pairing states, 
and half-fluxoid states. Some of these properties are intimately connected with the properties as a topo- 
logical superconductor The article concludes with a summary of knowledge emerged from the study of 
Sr2Ru04 that are now more widely applied to understand the physics of other unconventional supercon- 
ductors, as well as with a brief discussion of relatively unexplored but promising areas of ongoing and 
future studies of Sr2Ru04. 
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1. Spin-triplet superconductors 

1.1 Candidates of Spin-triplet superconductors 

In the last three decades, and particularly since the dis- 
covery of high-transition-temperature (high-Tc) superconduc- 
tivity of the cuprates, ' studies of "unconventional" supercon- 
ductivity have been one of the main topics in condensed- 
matter physics. Here we designate the term "unconventional" 
as the pairing based on non-phonon mechanisms.^ The un- 
conventional superconductivity is mainly found in heavy- 
fermion superconductors (since 1978),^ high-Tc cuprates 
(since 1986),' organic superconductors (since 1980),"* ruthen- 
ate superconductors (since 1994),^ and iron-pnictide super- 
conductors (since 2008).^ Unconventional superconductivity 
is characterized by the anisotropic gap function or order pa- 
rameter which is integrated to be zero or a small value due 
to the variations of the wave function "phase", in contrast to 
an ordinary i-wave state. In many of them, including high-r^ 
cuprates and iron pnictides, the electrons are clearly paired 
in spin-singlet states. In this point of view, they are similar 
to conventional i-wave superconductors, in which the spin- 
degrees of freedom is lost in the charged superfluids. Spin- 
triplet superfluid states are fully established in the Fermi liq- 
uid ^He,^ '^ for which spin and mass supercurrents emerge in 
the charge-neutral superfluids. The question is then whether 
or not spin-triplet supercoiiductors exist, and what novel su- 
perconducting properties they may exhibit due to their charge 
and spin supercurrents. 

There are several classes of candidates of spin-triplet 
superconductors represented in Table III. We should note 



here that due to strong spin-orbit interactions of heavy ele- 
ments, spin is not a good quantum number in the /-electron 
based heavy fermions; their possible pseudo-spin-triplet pair- 
ing may be better termed as "odd-parity" pairing for their or- 
bital wave function symmetry. Keeping this in mind, we will 
nevertheless use the term "spin-triplet" whenever it causes 
no confusion. Among the heavy fermions, UPt3 is the lead- 
ing candidate of a triplet superconductor^ The measurements 
of spin susceptibility indicate that pairing in UNi2Al3"^ is 
also triplet," whereas the pairing in UPdiAls is clearly spin- 
singlet.'^'^ Recent progress in the studies of ferromagnetic 
heavy-fermion superconductors such as UGe2,''* URhGe,'^ 
and Ulr,'^ and in particular UCoGe,'^ motivates researchers 
to seek for obtaining direct evidence for spin-triplet pairing. 

Among the organic superconductors, possible spin- 
triplet pairing in quasi-one-dimensional superconductors 
(TMTSF)2X (X = CIO4 and PFg) has been considered.^^-^' 
However, spin-singlet pairing is demonstrated by nuclear 
magnetic resonance (NMR) Knight shift for X = C104.^^ In 
addition, recent results of detailed field-orientation depen- 
dence of the onset temperature derived from resistivity in- 
dicate that the unusual survival of superconductivity in high 
magnetic fields is attributable to the formation of the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state of spin-singlet su- 
perconductivity.^'''' ^'^ 

A new class of spin-triplet pairing has been recognized 
in some superconductors having crystal structures without in- 
version symmetry, such as LiiPtsB.^"* It is believed that if the 
spin-orbit splitting of the Fermi surfaces due to internal elec- 
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Table I. Selection of candidates of spin-triplet superconductors. HF: heavy fermion superconductors, NCS: Noncentrosymmetric superconductors, FM: 
fen'omagnetic superconductors, *: superconductivity underpressure. 



Materials 


Classification 


Spin evidence of triplet pairing 


Properties 


^He 


Superfluid 


magnetization, NMR etc. ' 


p-wave, A phase is chiral 


Sr2Ru04 


Oxide 


NMR, polarized neutron 


2D analogue of~^He-A 








Chiral p-wave 


UPt3 


HF 


NMlT'*^ 


/-wave 


UBei3, URu2Si2, UNiaAlj 


HF 


NMRi^ 




UGe;, URhGe, UCoGe 


FM, HF 


Indirect 


Anomalous i/c2^^"^^ 


Ulr* 


NCS,FM,HF 


Indirect 




CelrSi; 


NCS, HF 


NMR23 




Li2Pt3B 


NCS 


NMR24 




CePtsSi 


NCS, HF 


Indirect 




CeRhSi; 


NCS, HF 


Indirect 


Anomalous i/c2^^ 


S/FM/S 


junctions 


Indirect (/c)^^-^^ 


Odd-freq., even-parity, i-wave 



trie fields exceeds the superconducting gap energy, the pairing 
state is no longer classified as purely spin-singlet or triplet, but 
their mixing occurs. •'^•■'^ As a result, a variety of unusual be- 
havior characteristic of spin-triplet pairing such as extremely 
high upper critical field //c2 emerges. Examples are CePtsSi,^' 
CeRhSis,^'^ etc. In particular, anisotropic spin susceptibility 
in the superconducting state has been reported for CelrSis.^^* 
As Mukuda et al. explained, if the spin susceptibility is domi- 
nated by strong Rashba-type spin-orbit coupling. Knight-shift 
measurements do not give direct evidence for the admixture 
between the Cooper-pairing states of even and odd parity. The 
anisotropy may well be dominated by that of Van-Vleck-like 
susceptibility. 

Sr2Ru04, the main subject of this review, is the first ox- 
ide superconductor in the same layered structure as high-Tc 
cuprates but with a low superconducting transition tempera- 
ture of 1.5 K.^' ^*^ The availability of high-quality single crys- 
tals and the relative simplicity of its fully-characterized Fermi 
surfaces promoted a large number of experimental as well as 
theoretical studies. Sr2Ru04 is now established as one of the 
archetypal unconventional superconductors. In particular, it 
is considered most probably as a spin-triplet superconductor, 
comparable with the odd-parity, pseudo-spin-triplet supercon- 
ductor UPtj."*" 

We should add that the emergence of spin-triplet pair- 
ing states is also discussed in superconductor-ferromagnet- 
superconductor (S/FM/S) junctions. By the proximity effect 
into ferromagnetic half-metal such as Cr02, i-wave pair am- 
plitude may penetrate as spin-triplet i-wave state. ^^'^^ Such 
an unusual state is known as odd-frequency pairing."" Re- 
cently, the long-range proximity effect due to odd-frequency 
spin-triplet currents is also demonstrated using ferromagnetic 
hybrid layers with some disorder and rotation of magnetiza- 
tion.28.29 



1.2 Unconventional superconducting properties of Sr2RuOi\ 
From its Fermi surface topography and anisotropy, 
Sr2Ru04 is characterized as a highly anisotropic supercon- 
ductor with the anisotropy ratio of 20. Nevertheless, the out- 
of-plane coherence length is several times greater than 
the interlayer distance to allow interlayer coherence. Ta- 
ble II summarizes basic superconducting parameters obtained 
mainly from the anisotropic //c2 and the specific heat. We note 
that the evaluation of assumes orbital limiting of H^2- 

Sr2Ru04 exhibits a number of phenomena reflecting its 
unconventional superconductivity. The extreme sensitivity to 
non-magnetic impurities and defects indicates non-i-wave 
pairing with sign-changes on the superconducting gap. In 
fact the transition temperature Tc diminishes at the impu- 
rity concentration at which the mean impurity-impurity dis- 
tance becomes comparable to the in-plane coherence length 
^o/,, consistent with the modified Abrikosov-Gorkov predic- 
tion. ''2''*'' Absence of coherence peak (Hebel-Schlichter peak) 
and the power-law temperature dependence in the NMR spin- 
lattice relaxation rate l/Ti indicate non-5-wave with highly 
anisotropic gap structure."''* The magnetic field dependence of 
the specific heat at low temperatures reveals a plateau in the 
intermediate H, indicative of the multi-gap superconductivity 
reflecting the existence of three Fermi surfaces with different 
orbital characters."'^ 

It should be mentioned that the theoretical proposal for 
the similarity between the superconductivity of Sr2Ru04 and 
the spin-triplet superfluidity of ''He soon after the discovery of 
ruthenate superconductivity"*^ strongly motivated the experi- 
mental efforts to determine the symmetry of unconventional 
superconductivity in Sr2Ru04. This leads to the first direct ex- 
perimental evidence of spin-triplet pairing in Sr2Ru04 by the 
measurements of the electron spin susceptibility with NMR."*^ 
Combined with the observation of internal magnetic field by 
(xSR,"*** this result yields a complete description of the super- 
conducting state of Sr2Ru04 as a two-dimensional analogue 
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Table III. Selection of candidates of spin-triplet superconductors. HF: heavy fermion superconductors, NCS: Noncentrosymmetric superconductors, FM: 
fen'omagnetic superconductors, *: superconductivity underpressure. 



Materials 


Classification 


Spin evidence of triplet pairing 


Properties 


^He 


Superfluid 


magnetization, NMR etc. ' 


p-wave, A phase is chiral 


Sr2Ru04 


Oxide 


NMR, polarized neutron 


2D analogue of^He-A 








Chiral p-wave 


UPt3 


HF 


NMlT'*^ 


/-wave 


UBei3, URu2Si2, UNiaAlj 


HF 


NMRi^ 




UGe;, URhGe, UCoGe 


FM, HF 


Indirect 


Anomalous i/c2^^"^^ 


Ulr* 


NCS,FM,HF 


Indirect 




CelrSi; 


NCS, HF 


NMR23 




Li2Pt3B 


NCS 


NMR24 




CePtsSi 


NCS, HF 


Indirect 




CeRhSi; 


NCS, HF 


Indirect 


Anomalous i/c2^^ 


S/FM/S 


junctions 


Indirect (/c)^^-^^ 


Odd-freq., even-parity, i-wave 



Table II. Superconducting parameters for Sr2Ru04. //c is the thermody- 
namical critical field, f is the coherence length, A is the penetration depth, 
and K is the Ginzburg-Landau parameter. 'Conventional orbital depairing 
for Hc2 is assumed. 



Parameter 






ab 


c 


Tc 


(K) 


1.50 








(T) 




1.50 


0.075 




(T) 


0.023 






m 


(A) 




660 


33* 


m 


(A) 




1900 


3.0 X 10^ 








2.6 


46 






20 







of the ^He-A phase, i.e., a spin-triplet chiral p-wave state. 

1.3 Scope of this review 

There have been several review articles devoted to the 
superconductivity of Sr2Ru04.^^'^"^'^'^'^' In addition, we men- 
tion a review on the Fermi-liquid properties in the normal 
state^^ and stories behind the discovery of its superconductiv- 
jjy 53,54 Especially, an extensive review of the developments 
up to the spring of 2002 by Mackenzie and Maeno has been 
cited more than 500 times. Since then, much research progress 
has been made on Sr2Ru04 in both experiments and theories, 
especially concerning the superconducting gap structure, the 
order-parameter symmetry, and the properties of supercon- 
ducting states injunctions and at boundaries. 

The purpose of this review is to make a critical exami- 
nation of the spin-triplet pairing scenario of Sr2Ru04, in or- 
der to fully understand the behavior of such superconducting 
state. The ideas emerged from the studies of Sr2Ru04 include 
multi-gap superconductivity of multi-band systems, physics 
of superconductivity with time-reversal symmetry breaking, 
pairing mechanism in highly correlated materials in which 



spin fluctuation is not necessarily strong. Such analysis should 
help understanding the properties of other unconventional su- 
perconductors, and moreover help resolving some of the long 
standing remaining issues of spin-triplet superfluid ^He.'*" 

In § 2, we first introduce the order parameter of spin- 
triplet superconductivity, the c/-vector, especially the c/-vector 
states relevant to the superconductivity of Sr2Ru04. We ex- 
amine the experimental evidence for spin-triplet pairing based 
on traditional probes, but with innovative approaches. A vari- 
ety of NMR Knight-shift experiments and polarized neutron 
experiments establish the spin-triplet pairing. We also address 
alternative possibilities of the li-vector orientation consistent 
with the NMR. We then describe experiments establishing the 
time-reversal symmetry breaking, or chiral nature, of the or- 
der parameter. 

Section 3 summarizes the theoretical advancement in un- 
conventional pairing mechanism by Coulomb repulsion. The 
importance of non-boson exchange, namely vertex-correction 
processes, is emphasized. We also review and examine other 
approaches to the unconventional pairing mechanism for 
Sr2Ru04. The important roles played by the spin-orbit inter- 
action, which determines the direction of the d- vector, are also 
discussed. 

Section 4 deals with the superconducting properties of 
Sr2Ru04 which remain unresolved within the current theoret- 
ical models. These are the origin of the strong suppression of 
//c2 with the fields parallel to the Ru02 plane and the emer- 
gence of a second superconducting phase at low temperatures 
under the same condition. 

Section 5 describes novel superconducting phenomena in 
Sr2Ru04, such as enhanced superconductivity (the so-called 
3-K phase) at the interface between Sr2Ru04 and Ru in the 
eutectic system, and novel proximity effects in various su- 
perconducting junctions involving Sr2Ru04. Descriptions of 
Andreev bound states in terms of odd-frequency pairing and 
chiral edge states are introduced. Many of these phenomena 
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are characteristics of a topological superconductor We also 
describe the current experimental status of observing half- 
fluxoid states in micron-scale rings of Sr2Ru04 crystals and 
discuss the physics of half-quantum vortices in spin-triplet su- 
perconductors. 

Finally in § 6, we summarize the present status towards 
establishing the spin-triplet superconductivity in Sr2Ru04 and 
discuss the future prospects of the physics of topological su- 
perconductivity. 

2. Evidence for spin triplet pairing 

2.1 Order parameter of spin-triplet superconductivity: the 
d-vector 

We start with a useful description of the order param- 
eter of spin-triplet superconductivity, namely the tZ-vector. 
The superconducting order parameter is given by a function 
of spin and momentum in general: Aa-a'(k). Roughly speak- 
ing, A(j-o-'(^) is regarded as a wave function of a Cooper pair 
formed by two quasi-particles whose momenta and spins are 
{k, cr) and {—k, cr'). The magnitude of the total spin of a 
Cooper pair is 5 - 1 in spin-triplet pairing, in contrast to 
5 = in spin-singlet pairing. In spin-triplet pairing, the de- 
grees of freedom about the direction of the spin S = 1 re- 
main even in the superconducting state. For spin-triplet pair- 
ing, Aa-a-'{k) is expressed using the li-vector as:^^ 

K^,{k)^m{k)-(T)(r,:\^^,, (1) 

or 

/ A„(^) ^^^{k) \ 

\ ^y^{k) Au(A:) ] 

-dAk) + id,{k) d,{k) \ 
\ d,(k) dAk) + idy{k) j' 

where cr is the Pauli spin matrix vector: (r - (ctv, cr,,, cr.). 
The order parameter of triplet pairing states is described com- 
pletely by the three components of the c/-vector d(k). d(k) is 
an odd-parity function of k. The c/-vector behaves as a vector 
under rotations in spin space. Here we should note that the di- 
rection of the c/-vector is perpendicular to the total spin 5 of a 
Cooper pair. To see this simply, we consider the case of d\\z, 
i.e., d(k) = (0, 0, c/;(A;)). Then we have An(A;) = A^Cit) = 0, 
i.e., S- - 0, which means d ± S. In general, each compo- 
nent of the d-vector {dx, d^, or d~) represents the pair am- 
plitude with the Cooper-pair spin perpendicular to the cor- 
responding basis (x, y, or z). Triplet pairing states satisfying 
d(k) X d*(k) - (d(k) x d*(k) + 0) are called unitary (non- 
unitary) states. In non-unitary states, spins of the Cooper pairs 
are polarized. We can see this simply by considering the quan- 
tity |An(A;)p - |Ai|(jt)|2 = T\{d{k) x + for non- 
unitary states. 

The order parameters of superfluid phases of ^He are ex- 
pressed as:^^ d{k) oc z{kx±\ky?) (A phase), d{U) oc xk x+yk.,.-^zk- 
(B phase), d{k) cc (x + \y)(kx + i^y) (Al phase emerging un- 
der magnetic fields). Both of the A and B phases are unitary 
states, while the Al phase is a non-unitary state with spin po- 



larization. For the state d(]i) oc z{kx + i^y), we can choose a 
new spin coordinate system in which A|x(^:) = t^Y\{k) = 0, 
i.e., dAk) = at all k. To see this, we rotate the spin coor- 
dinate system by 2;7r/3 around the [111] axis. Then the basis 
vectors x, y and z are transfered to z, x and y, respectively, and 
consequently we have dAk) - d^ik) - for all k in the new 
spin coordinate system. In general, if we can choose appro- 
priate spin coordinates by which d,{k) = for all k, then we 
call the pairing state "equal-spin pairing" (ESP) state. The A 
and Al phases are ESP states, while the B phase is not. For 
those not familiar with the description of the spin-triplet or- 
der parameter, the review by Mackenzie and Maeno (§ IV and 
Appendix D), as well as the Appendix of this review, gives 
an introductory but more detailed treatment in terms of the 
c/-vector. For example, one can find the relation between the 
spin bases and the c/-vector bases, and the expression of the 
superconducting gap in terms of the d-vectoi. 

Down to the low temperatures, Sr2Ru04 maintains a 
tetragonal structure with the crystal point group symmetry 
Z)4h.^^ Neglecting the dispersion along the out-of-plane c di- 
rection, possible spin-triplet states are limited to those for 
the two-dimensional square lattice with C4v symmetry and 
are given in Table IV. Figures 1(a) and 1(b) depict the Fi 
and Fs states, respectively. The former can be thought of as 
the quantum-mechanical superposition of the two states with 
antiparallel combinations of spin and orbital wave functions, 
with quenched overall orbital angular momentum: 

d(k) oc (x -H iy)(kx - iky) -i- (x - iy)(kx -H i^,.) (3) 

The li-vector x + iy lies in the basal plane, which for illus- 
trative purpose is depicted by yellow arrows in Fig. 1(a). The 
Fj states in Table IV are called the chiral pairing states, be- 
cause these states possess two polarizations of relative orbital 
angular momentum of pairing quasi-particles: left- and right- 
handed polarizations correspond to kx+\ky and k^-iky, respec- 
tively. They are the states with the orbital angular momentum 
= -1-1 or -1 as represented by big arrows in Fig. 1(b); the 
Cooper-pair spins he in the basal plane asd ± S. We note that 
all five states Fi - F5 are the so-called ESP states, for which 
the spin state of the Cooper pair is always a linear superpo- 
sition of S - - +1 and 5 , = - 1 states by a suitable choice of 
spin axes. 

2.2 Spin part of the order parameter 

The direct evidence for spin-triplet pairing in Sr2Ru04 
is based mainly on the following two kinds of experiments. 
One is the electronic spin susceptibility measurements by the 
NMR Knight shifts of both '^O and ''''Ru nuclei.^^ ^s The 
other evidence comes from the measurements of local magne- 
tization by polarized neutron diffraction; the so-called ShuU- 
Wedgewood experiments.^"^ 

In this subsection, we focus on the NMR experiments in 
some detail. The NMR Knight shift measures the effective 
field at the nucleus produced by the electrons, and thus is re- 
lated to the microscopic spin susceptibility at the NMR nu- 
clear site. Since the spin susceptibility in the superconducting 
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Table IV. List of possible five rf-vector states for Sr2Ru04. Irreducible representations of the C4,. symmetry group are displayed by both of the Bethe and 
Mulliken notations in the first two columns. The minus sign in the Bethe notation and the index "u" (denoting ungerade) in the Mulliken notation represent 
odd parity. TRSB in the fifth column means time-reversal symmetry breaking. L^ and Fj are analogous states to the ^He-B and ^^He-A phases, respectively. 
Lj corresponds to the chiral p-wave state. In the absence of the spin-orbit interaction among Ru-4rf electrons, i.e., assuming isotropy in spin space, any 
superpositions of these five states are also possible. Note that the column "rf-vector" displays only the symmetry of d(k) as a function of k. For example, 
if we take the functions sin/:, and sin/:,, for the i,-- and iry-symmetry, then we have d{k) oc zi^mk^ ± i sin/:,,) for Fj . All of the five states are unitary ESP 
states, and give full energy gap at the Fermi level. 



Bethe 


Mulliken 


(i-vector 


direction of d 


TRSB 


Analogy to ^He 




Au, 


xk^ + yky 


d II ab 


No 


BW state (B phase) 




Aiu 




d II ab 


No 






B\u 


xk^ - yky 


d II ab 


No 






Biu 


xky + ykjc 


d II ab 


No 






Eu 


z(kj, + iky) 


d\\c 


Yes 


ABM state (A phase) 



state cannot be extracted from the bulk susceptibility due to 
the Meissner shielding, the Knight-shift measurement is con- 
sidered as the most reliable method to measure the spin sus- 
ceptibility in the superconducting state.'*"* 

The Cooper pair spin susceptibility of Sr2Ru04 can be 
reliably extracted from the NMR Knight shift K for the fol- 
lowing reasons. First, the NMR and nuclear quadrupole reso- 
nance (NQR) measurements have been performed on a vari- 
ety of nuclei, three different crystal sites of '^O (NMR and 
NQR), '''^Ru (NMR), and '^'Ru (NQR), covering wide pa- 
rameter range in the field-temperature (H-T) phase diagram. 
Thus accidental signal cancellation due to the structure factor 
is avoided. In addition, it is noticeable that one can also ob- 
tain the bulk susceptibility and nuclear spin-lattice relaxation 
rate (l/Ti) within the same measurement setup. This fact en- 
ables us to carefully check the consistency of the data from 
different viewpoints. Second, overheating of the sample by 
RF pulses is avoided by monitoring the intensity of the NMR 
peak itself, thereby using the NMR nuclei as an internal ther- 
mometer (Supplement of Ref. 47). Third, based on the sys- 
tematic studies of K among different ruthenates, i.e., RUO2, 
Sr2Ru04, and CaRuO^, it is clear that a large negative value 
of K for Sr2Ru04 originates from the contribution from core 
polarization due to li-electrons (Fig. 2).^" From such analy- 
sis, the contribution from the orbital part is estimated as 1.0% 
while the spin contribution is -4.5% for A" || c of Sr2Ru04. 
With the anisotropy of the spin part of only 2% as extracted 
from thtXab ^^Xc plot,^' the dominance of the spin part in the 
susceptibility^^' ^-^ is illustrated in Fig. 3. Such properties of 
Sr2Ru04 ensure that its NMR Knight shift serves as a reliable 
high-sensitivity probe for the Cooper pair spin susceptibility. 

Ishida and coworkers have measured Knight shifts K{T) 
at the Ru and O sites in high-quality single crystals. Mea- 
surements were performed using a ^He-'^He dilution refrig- 
erator with the sample crystals directly immersed in liquid 
helium to ensure good thermal contact to the bath. As for 
the '^O Knight-shift measurements, they assigned NMR sig- 
nals arising from the different crystal sites (planer and api- 
cal O sites), and estimated the spin part of the Knight shift 
above at each site. They investigated temperature depen- 



dence of the Knight shift at the planer O site, where the spin 
density is 6 times larger than that in the apical O site, and 
showed that its Knight shift remains unchanged in the field 
range (0.35 T < hqH < 1.1 T) applied parallel to the RUO2 
plane.'*' 

They also measured the Knight shift at the Ru site in 
the RUO2 plane. The advantage of the Ru Knight shift mea- 
surements is that the value of the ^^Ru Knight shift is an 
order of magnitude larger than that of the '^O Knight shift, 
since the ^''Ru shift is strongly affected by Ru-4c/ electronic 
spins through the large hyperfine coupling constant of Acp — 
-25 T/yt/B originating from the core polarization effect, where 
yt/B is the Bohr magneton. This means that the Ru Knight shift 
can detect much smaller changes of the Ru-4^f electronic state 
when superconductivity sets in. The temperature dependence 
of the ^''Ru shift was measured in the fields of 0.68, 0.9, and 
1.05 T parallel to the RUO2 plane and the spin part of the 
Knight shift C'^K^^i^ = K{T)-Ko,b) obtained at 0.9 T is shown 
in Fig. 4. 

^^^spin does not change on passing through Tc ~ 1 .0 K at 
/ioH ~ 0.9 T. If a spin-singlet d-wave state with a line-node 
gap were realized, the T dependence of "^^A^spin would behave 
as drawn by the dashed curve in Fig. 4. Here we adopted the 
magnitude of the gap 2Ao = 4^b7^c, which has been used for 
analyses of various physical quantities. It is obvious that the 
unchanged Knight shift in the superconducting state at both 
the O and Ru sites in the RUO2 plane cannot be understood by 
the spin-singlet scenarios, but strongly suggests the Cooper- 
pair spins are in the triplet state and can be polarized to the 
applied field parallel to the RUO2 plane. 

In general, the spin susceptibility in the spin-triplet su- 
perconducting state is expected to show anisotropic temper- 
ature dependence reflecting the magnitude and direction of 
the applied field. When the interaction acting on the triplet 
pairs is strong enough to lock the spins of a superconduct- 
ing pair in a certain direction of the crystal, the anisotropy 
of the spin susceptibility, i.e., the anisotropic behavior of the 
Knight shift should be observed in the superconducting state. 
The Knight shift is unchanged when magnetic fields are par- 
allel to the spin direction (d ± H). But when magnetic fields 
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Fig. 1. (Color online) Illustrations of Cooper-pair spins and orbital angular 
momenta for spin-triplet superconductivity. The Cooper-pair spins S are 
depicted by small arrows in various colors; the rf-vector shown by yellow 
arrows is defined perpendicular to the spins. The orbital angular momen- 
tum L of a Cooper pair is shown by a large blue arrow, (a) The S and 
L vectors for the Fi state, d(k) = Ao(xkji + yky). It can be thought of as 
a quantum-mechanical superposition of Sy = +1 and Sj = -1 states (see 
text), (b) The promising ground state d{k) = z&o(kt +iky) with the rf-vector 
pointing in the c direction, (c) Possible states under magnetic fields. The 
c(-vector rotates within the basal plane in response to the external magnetic 
fields. The Cooper-pair spins can be polarized by magnetic fields of any 
directions. 



are applied perpendicular to the spin direction {d \\ H), the 
Knight shift decreases depending on the superconducting gap 
structure. On the other hand, if the energy of the pinning in- 
teraction is weaker so that the spins can reorient easily toward 
the applied magnetic field, the Knight shift is invariant in any 
field direction. Thus the Knight shift measurements in small 
fields are important to understand the pinning interaction of 
the c/- vector. 

In order to measure the Knight shift in such small fields, 
Murakawa et al. employed a ""Ru NQR spectrum for ob- 
serving the signals, and measured the Knight shift at the Ru 
gjjg63,64 shown in Fig. 5(a). They measured the tempera- 
ture dependence of the resonance frequency of the Ru NQR 
signal affected by the applied field parallel to the Ru02 plane. 
The expected decrease of the spin susceptibility in the case of 
the spin-singlet pairing or spin-triplet with the c/-vector in the 
RuOa plane is shown by the dotted curve. Obviously such sup- 
pression was not observed, and the spin susceptibility is in- 



Fig. 2. Knight shift vs bulk susceptibility determined in the range T = 1.4 
- 20 K for a weakly-correlated metal Ru02, an electron-con'elated super- 
conductor Sr2Ru04, and a correlated, nearly ferromagnetic metal CaRuOs . 
The hyperfine coupling constant due to the inner core polarization A^p ~ 
-22.2 T//iB is estimated from a lineal' slope with the assumption that 
Korb = 1.59% is common among these ruthenates. The value of A^p is 
close to Ahf ~ 30 T//jb in the FM state of SrRuOs . A large negative K 
value of Sr2Ru04 indicates that its susceptibility is dominated by the d- 
electron spin susceptibility (Fig. 5 in Ref. 60). 
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Fig. 3. (Color online) Bulk susceptibility of Sr2Ru04 (Refs. 62 and 63) 
illustrating the spin and orbital contributions (see text for details). 



variant through instead. The same measurement was done 
in field as small as 55 mT confirming the invariance of the 
spin susceptibility. The occurrence of the superconductivity 
in the field was confirmed by the measurements of the Meiss- 
ner signal as well as l/Ti of Ru using the identical set up of 
the Knight-shift measurements, as shown in Fig. 5(b). Note 
that the Knight-shift spectra were obtained within 100 |.is af- 
ter the RF pulse with the total duration of 100 [is. In contrast, 
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Fig. 4. Knight sliift of "Ru nuclei of Sr2Ru04 with the applied field par- 
allel to the basal ab plane (Ref. 62). The spin part K,{T), after subtracting 
the temperature-independent orbital part K^^f, from the total shift K(T), is 
plotted. 
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Fig. 6. ""Ru Knight shift ""/sT,. of Sr2Ru04 (a) at 90 mT for H \\ c as 
small as 20 mT and (b) at various temperatures at 44 mT (Ref. 64). One 
can see in (a) that the orbital contribution to is -1.6%. In (b) the data 
point in a triangle is obtained by intentionally varying the external field by 
1 mT. 



I IT I was obtained within 100 ms after the pulse. Absence 
of excessive heating by the RF pulses for the NMR measure- 
ments was carefully confirmed by examining the temperature 
dependence of the spectral peak intensities. 

The Knight-shift measurements in magnetic fields per- 
pendicular to the RuOa plane (parallel to the c axis), denoted 
as Kc, were performed in the same way. However the mea- 
surements in this direction is difficult for the following two 
reasons: (1) H^i along the c axis is only 75 mT, 20 times 
smaller than //c2 along the RuOa plane; (2) the supercon- 
ducting fraction with penetrating magnetic field is much more 
suppressed since the Ginzburg-Landau parameter is as small 
as 2.6 for H \\ c. Thus, the measurements for H \\ c were 
performed with utmost care. 

Figure 6(a) shows Kc measured at T =90 mK in the 
field range from 20 mT to 120 mT, and Fig. 6(b) shows the 
temperature dependence of Kc measured at 44 mT. In both 
measurements, no appreciable changes can be detected be- 
yond the experimental error The solid curves in Figs. 6(a) and 
6(b) are conceivable dependences, when the superconducting 
c/-vector is locked along the c axis. Such dependences are in- 
consistent with experimental results: The spin susceptibility 
is unchanged also for magnetic fields parallel to the c axis. 

Figure 7 summarizes the temperature and field ranges of 
the Knight-shift measurements performed so far; no apprecia- 
ble suppression has been detected in any of the measurements. 
One of the plausible interpretations is that the c/- vector orients 
perpendicular to the external fields greater than 20 mT, irre- 
spective of the field direction. 



Fig. 5. (Color online) (a) Low-field Knight shift of Ru nuclei of 
Sr2Ru04 with the field parallel to the basal ab plane. The Knight shift 
was obtained from the specfi'al shift of the NQR frequencies in the pres- 
ence of a magnetic field, (b) Ac susceptibility Xin^ and nuclear spin-lattice 
relaxation rate l/Tj obtained in the identical setup as used for Fig. 5(a). 
Taken from Ref. 63. 



2.3 Chiral order parameter 

The odd parity of the orbital part of the order pa- 
rameter has been unambiguously demonstrated by phase- 
sensitive measurements by Nelson et al.^^-^^ The critical cur- 
rent of Auo.5lno.5-Sr2Ru04 superconducting quantum inter- 
ference devices (SQUlDs) containing a loop of conventional 
s-wave superconductor Auo.slno.s and Sr2Ru04 prepared on 
Sr2Ru04 single crystals was found to be a maximum for de- 
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Fig. 7. (Color online) Field-temperature phase diagram of Sr2Ru04 sum- 
marizing the temperature or field sweeps of the Knight-shift measurements 
reported so far (Ref. 63). 



vices with junctions on the same side of the crystal and a min- 
imum for devices with junctions on opposite sides, in the limit 
of zero magnetic flux as shown in Fig. 8. These findings indi- 
cate that the phase of the superconducting order parameter in 
Sr2Ru04 changes by n under inversion. This result verifies the 
odd-parity pairing symmetry corresponding to the formation 
of spin-triplet Cooper pairs in Sr2Ru04. 
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Fig. 8. A ff-junction SQUID consisting of Auo.sliio.s and Sr2Ru04.''' The 
minimum of its critical current extrapolates to zero towards the supercon- 
ducting of an .s-wave superconductor Auo.slno s, where the induced flux 
due to the asymmetry of the SQUID is minimized. 



Among the candidate odd-parity spin-triplet pairing 
states listed in Table IV, the Fg state with the orbital order 
parameter A(^) - Ao(A:.v + iky) is the so-called "chiral super- 
conducting state", which is accompanied by the broken time- 
reversal symmetry (TRS). The term "chiral" is used for the 
TRS broken state associated with the orbital part of the wave 
function. The TRS breaking was first demonstrated by the 
emergence of internal magnetic fields observed in the muon 
spin resonance experiments.'*** Since the specific heat demon- 
strates no residual density of states (DOS), the broken TRS 
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Fig. 9. Rotation angle A6k of the reflected light from the surface of a 
Sr2Ru04 crystal due to magneto-optic Kerr effect (MOKE). The change 
of about 60 nrad below is attributed to the broken time-reversal sym- 
metry of chiral superconductivity. The sign of A0k can be manipulated by 
the direction of the external field applied before the measurements, inter- 
preted as the alignment of chiral superconducting domains. From Xia et 
al. (2008).*** 



State associated with the spin part, the so-called "non-unitary 
state", is unlikely realized in Sr2Ru04 at least under zero mag- 
netic field. 

An important progress has been made by the devel- 
opment of ultra-high sensitive magneto-optic Kerr effect 
(MOKE) apparatus based on a Sagnac interferometer^^ The 
rotation of the polarization plane of the reflected light below 
Tc shown in Fig. 9 indicates the broken TRS in the super- 
conducting Sr2Ru04.^^ It is shown that the Kerr signal does 
not couple strongly with the vortices and that the sign of the 
Kerr signal is reversed by the application of the field over 
about 5 mT before the measurements under zero-field are per- 
formed. The latter is attributed to the need for training the 
chiral domains, i.e., enlarging the size of the chiral domain 
having the same sign of the orbital angular momentum of the 
Cooper pairs as the external field. Considering the size of the 
focused laser light used, the trained chiral domain size is es- 
timated to be on the order of a few micrometers. It is worth 
mentioning that the new technique of MOKE experiments has 
successfully been applied to investigate broken TRS associ- 
ated with the inverse proximity effect in the bilayer films con- 
sisting of an i-wave superconductor and a ferromagnet.^'' 

The MOKE measurements for Sr2Ru04 stimulated a 
number of theoretical investigations. ^"^"^^ In general, the Kerr 
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rotation angle is expressed in terms of the ac Hall conductivity 
ax\'(cL>), i.e., the off-diagonal part of the ac conductivity tensor. 
Using appropriate parameters for Sr2Ru04, Yakovenko eval- 
uated the Kerr rotation angle A6k for the chiral p-wave pair- 
ing state, and obtained A6k ~ 230 mad.™ Goryo presented a 
mechanism of impurity-induced Kerr rotation. He calculated 
the Hall conductivity within the low-order expansion in the 
impurity strength, and found that the leading contribution is 
given by similar diagrams to the skew-scattering diagrams in 
the extrinsic anomalous Hall effect. In this impurity-induced 
mechanism, the MOKE would be suppressed or would be 
zero for any unitary pairing states other than the chiral pairing 
state. The rotation angle A6k is evaluated to be about 30 nrad 
for a realistic impurity mean distance i = 5000 A, in agree- 
ment with the experimental values.^' These studies on MOKE 
suggest that the chiral pairing state is realized in Sr2Ru04. 

The chiral superconducting state of Sr2Ru04 is supported 
by a number of experiments, including yuSR, Kerr effect, and 
Josephson effect. The chirality, or the TRS breaking in the or- 
bital wave function, is believed to originate from the kjc + iky 
wave function associated with the orbital angular momentum 
L~ = 1 of Cooper pairs. In the bulk of Sr2Ru04 crystals, it is 
expected that domain structures of superconducting regions 
with opposite chirality form. Strong evidence for the pres- 
ence of such order parameter domains as well as for the do- 
main wall motion were obtained through anomalous behavior 
of diffraction patterns of Josephson junctions on single faces 
of the crystals by Kidwingira et aV^ The observed telegraph 
noise in the critical current as a function of magnetic field or 
time, responsible also for hysteresis observed in field sweeps 
of the critical current, is attributable to transitions between the 
chiral states of a domain or the motion of domain walls sep- 
arating them. Examining the modulation envelopes the data, 
they estimated average domain width of ~ 1 \\n\. The presence 
of such domains confirms the broken time-reversal symmetry 
nature of the superconducting pairing state in Sr2Ru04. 

At a boundary of a chiral superconducting mono-domain, 
net superconducting edge current is expected within a healing 
length proportional to the superconducting coherence length. 
However, this edge current is expected to be screened by the 
Meissner current extending to the penetration depth from the 
edge. Since the Ginzburg-Landau parameter is k - 2.6, this 
cancellation must be rather significant. Nevertheless, within 
the parameters of Sr2Ru04, it is expected that the net mag- 
netic field due to chiral edge current is as large as 1 mT. On 
the cleaved surface of a single crystal of Sr2Ru04, it is ex- 
pected that the chiral edge current exists at the bulk edge as 
well as along the boundaries of + iky (L- = +1) and k^ - iky 
(Lj. = -1) domains. 

Aiming at the direct observation of the magnetic field due 
to the chiral edge current, measurements of scanning SQUID 
microscope were performed. ^^'^^ In both experiments, spon- 
taneous magnetic field was not observed at the bulk bound- 
aries or as an internal pattern suggesting the presence of chi- 
ral domain boundaries. After various domain configurations 
are considered, it is concluded that for the domain size of 2 
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microns or greater, the magnitude of the magnetic field has to 
be less than 1 % of that expected. Another approach in search 
for the chiral edge current was taken by a torque magnetome- 
try on a micron-size Sr2Ru04 ring; but the result indicated no 
edge current of the expected magnitude. 

We note that Dolocan et al. also made direct imaging of 
magnetic flux structures using a scanning micro-SQUID force 
microscope, anticipating preferential penetration as well as 
pinning of vortices along the domain walls. They observed 
that as field parallel to the c-axis is increased, individual vor- 
tices coalesce and form flux domains. These observations im- 
ply the existence of a mechanism for bringing vortices to- 
gether overcoming the conventional repulsive vortex-vortex 
interaction. Simulation study of the magnetization process 
was performed by Ichioka et. al. for the multidomain state 
in a chiral p-wave superconductor, using the time-dependent 
Ginzburg-Landau theory.^^ They derived a vortex sheet struc- 
ture forming along a domain wall, and at higher fields the mo- 
tion of the domain walls so that the unstable domains shrink 
to vanish. Therefore, the single domain structure would be re- 
alized at higher fields. 

To explain the absence of magnetic field associated with 
the edge current, some scenarios have been proposed. Raghu 
et al. proposed an alternative triplet pairing state, '^'^ focus- 
ing on the hidden quasi-one-dimensional character of the xz 
and yz bands and the interplay between spin and charge fluc- 
tuations on these bands. They obtained a superconducting 
state without robust chiral edge modes along the boundary, 
thereby explaining the absence of experimentally detectable 
edge currents. Based on the Coulomb interaction mechanism 
presented in § 3, Tada et al. showed that the TRS can be re- 
stored near the (001) surface because of a Rashba-type spin- 
orbit interaction.^' They suggested that such pairing state at 
the interface is a promising candidate for the recently pro- 
posed time-reversal invariant topological superconductivity 
carrying helical edge current. The surface states they pro- 
pose naturally lead to cancelled charge current without mag- 
netic field around the edge, although the authors did not ex- 
plicitly mention a relation to the above experimental results. 

These null results for edge currents may also be related 
to a long-standing profound question concerning the size of 
the "intrinsic angular momentum" of chiral p-wave super- 
fluid. The chiral edge current originates from a topologi- 
cal property of the bulk and is related to the angular mo- 
mentum of Cooper pairs. ^'^ In a superconductor, Meissner 
screening current is spontaneously generated in response to 
the topological edge current. For ^He-A phase with A(^) - 
Aa(A:v + iky), the expected topological edge mass current has 
never been observed. The magnitude of the mass current de- 
pends on the magnitude of the total angular momentum (L) - 
{Nh/2){A/Ep)" . Depending on whether all the fermions con- 
tribute their pair angular momenta or only those in the interac- 
tion shell contribute them, the value of a varies as a = 0, 1 or 
2. As Leggett notes, this problem "is more than 30 years old 
and still has apparently not attained a universally agreed reso- 
lution."*" The estimated magnitude of the chiral edge charge 
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current in Sr2Ru04 has been largely based on the assumption 
of a = 0, but a reduction factor with (A/Ep) ~ 10^^ may have 
to be included. 

2.4 Superconducting gap structure 

At the time when the previous review by Mackenzie and 
Maeno was written, the main issue concerning the gap struc- 
ture was apparent controversy between the observed low- 
lying quasi-particle excitations and the predicted fully-gapped 
state. The former was deduced from the T-linear electronic 
specific heat coefficient^^ C/T ~ T and the T-cubed nuclear- 
spin relaxation rate^^ l/Ti ~ T^, as well as from the ther- 
mal conductivity, '^^^ penetration-depth,*^^ and ultrasound at- 
tenuation. The anticipation of the full gap is a consequence 
of the spin-triplet pairing state on the quasi-two dimensional 
Fermi surface in the tetragonal symmetry, d{k) = zAo(kt+iky) 
with the isotropic full gap \A{k^, = Aq. It was then rec- 
ognized that multiband effects are important: a large super- 
conducting gap corresponding to opens up on the "ac- 
tive" y Fermi surface derived mainly from the dxy orbitals of 
the Ru-4ii electrons and having 57% of the total electronic 
DOS, while the "passive" a and /? Fermi surfaces derived 
mainly from the dyz and d-x orbitals host a gap substantially 
smaller compared with the energy scale of Tc- The plateau 
in C/r with increasing H^ab is ascribable to the presence of 
the small gap in the passive bands with the DOS amounting 
to 43%. Moreover, nearly-zero gap is expected for these pas- 
sive bands. At the same time, various models have also been 
proposed which contain "horizontal" line nodes, with the gap 
zero at a certain k.. 

A solution to these apparent controversies has been 
largely given by the specific heat under magnetic fields ori- 
entated accurately to the crystalline axes.^^ '" Before describ- 
ing details of the superconducting gap structure investigation 
from these field-angle controlled specific-heat measurements, 
we first review the theoretically predicted orbital-dependent 
gap structure. At present, the gap structure semi-quantitatively 
consistent with the existing experiments is the one predicted 
by Nomura and Yamada.^^ As shown in Fig. 19, it consists of 
a large and anisotropic full-gap on the y Fermi surface, hav- 
ing a gap minimum along the F-M direction (corresponding to 
the Ru-O [100] direction), and of smaller gaps on the a and /3 
Fermi surfaces, amounting to 1/3 - 1/5 of the gap size on the y 
Fermi surface, and having a strong gap minimum accidentally 
almost zero along the F-X (diagonal) directions. 

The origins of such "orbital-dependent" gap structure are 
attributable to the differences in orbital symmetry of each 
band and in the resulting contribution to the pairing mecha- 
nism of each band. Such multiband gap structures with dif- 
ference in size and anisotropy, introduced in connection with 
Sr2Ru04,^"^ are proven to serve as an important concept also 
in the superconductivity of other systems such as MgEa**^ and 
iron-pnictides.''^ 

We next describe somewhat in more detail how the 
anisotropic gap structure has been deduced from the specific 
heat."*^ The quasi-particle spectroscopy using the specific heat 




-A A E-Ep 




Fig. 10. (a) Volovik effect: the energy shift of a quasi-particle spectrum 
due to the Doppler effect. Quasi-particle excitations around the nodes or 
minima of the superconducting gap for the applied magnetic fields parallel 
to (b) the nodal (or gap minimum) direction and (c) anti-nodal (or gap 
maximum) direction. 



under oriented magnetic fields bases its principle on the so- 
called Volovik effect.^* In the mixed state, the supercurrent 
velocities around the vortices are perpendicular to the applied 
field. Since the quasi-particle excitation spectrum is deter- 
mined with respect to the rest frame of the condensates, the 
presence of the supercurrents results in the Doppler-shift in 
the quasi-particle spectrum in the laboratory frame as illus- 
trated in Fig. 10(a): E{k) - Eo{k) + nif-Vpik) ■ Vs, where E{k) 
and £'o(^) are the quasi-particle spectra in the frames of the 
laboratory and the condensate, ;7Ze is the electron mass, Vf and 
Vs are the Fermi velocity and the supercurrent velocity, and 
k is the quasi-particle wave vector in the laboratory frame. 
For quasi-two-dimensional superconductors, Vf is expected 
to be parallel to k. As illustrated in Fig. 10(b), this Doppler 
shift creates field-induced quasi-particle excitations in the re- 
gions of the A:-vector around gap nodes or gap minima. Such 
field-induced quasi-particle excitations do not occur for the k- 
vector parallel to the applied field since Vf is perpendicular to 
Vs. For systems with four line nodes on the two-dimensional 
Fermi surface, for example, the total quasi-particle excitations 
are minimized for the field direction along the nodes. Since 
the specific heat divided by T or thermal conductivity divided 
by T is proportional to the quasi-particle DOS at low temper- 
atures, it takes a minimum if the field points to the direction 
of a nodal or minimum gap position. 

An apparatus developed to investigate the field-angle 
quasi-particle spectra is depicted in Fig. 11.^^ In this appara- 
tus, the so-called "vector magnet" consisting of split-pair su- 
perconducting magnet providing the horizontal field and the 
solenoid magnet providing the vertical field is placed in the 
dewar. This dewar rotates horizontally on the platform while 
the dilution refrigerator inserted in the dewar is fixed. This ap- 
paratus allows field alignment with respect to the crystalline 
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Fig. II. Apparatus used to measure anisotropic superconducting properties 
under magnetic fields of accurately and precisely aligned directions. From 
Deguchi et al. (2004).'"* 



axes of the sample with high accuracy and precision: the pre- 
cision in the azimuthal angle of ~ 0.001° and that of the 
polar angle of AG ~ 0.01° are readily achieved. Because we 
do not need to rotate the refrigerator, the thermal stability of 
the measurements is greatly enhanced. 

Figure 12 represents C/T with the magnetic field rotating 
in the basal plane. At relatively high temperatures and near 
//c2, the four-fold oscillatory component of C/T with mini- 
mum along the [110] direction is governed by the anisotropy 
of Hc2 with the maximum in that direction below 0.8 K.^^- 
The data at 0.31 K and 0.12 K clearly exhibit a reversal of 
the C/T oscillation when the field is reduced away from Hc2- 
In the Doppler shift analysis, the minimum along the [100] 
direction in this low-// and low-T region indicates that the 
minimum in the main superconducting gap is in the [100] di- 
rection. 

Curiously, as the field is reduced below a certain value 
at the lowest measurement temperature, the oscillatory ampli- 
tude in C/T diminishes (Fig. 13), contrary to the expectation 
of growing relative amplitude for a single-band case. Since 
the threshold field of diminishing oscillation corresponds well 
with the exhaustion of the quasi-particle excitations corre- 
sponding to the passive a and /3 bands, it is most naturally 
explained by the competition of the oscillations between the 
active and passive bands with opposite gap anisotropics. The 
specific-heat data under fields tilted from the basal plane as- 
sures this interpretation.^' 

In addition to the anisotropics due to H^i and nodal gap 
structure, recent experimental and theoretical studies revealed 
that the in-plane anisotropy of the specific heat or thermal 
conductivity in nodal superconductors can change sign, de- 
pending on the temperature and field range, due to the vortex 
scattering effect under rotating magnetic field. In fact, such 
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Fig. 12. (Color online) (left) Fourfold oscillation of specific heat of 
Sr2Ru04 under field rotation within the basal plane. Oscillations at high 
fields are attributable to the anisotropy of //c2-" In the low-temperature 
and low-field region, the oscillation emerges with the opposite phase from 
the Hc2 anisotropy. (right) Superconducting gap anisotropy on the active y 
band, deduced from the specific-heat oscillations. The gap minima are in 
the [100] directions. From Deguchi et al. (2004)."^ 



a sign change was observed at J ^ 0.1 Tc in C/T of the d^i^yi- 
wave superconductor CeCoIns.'"^ For Sr2Ru04, the observed 
large anisotropy at low temperatures probably represents the 
gap anisotropy of the main band since the measurements were 
extended to below 0.1 Tc, and the reduction of the anisotropy 
at the lowest temperature of the measurements is attributed to 
the multiband effect. 

The gap anisotropy consistent with these observations 
has been predicted by the microscopic theory based on the 
three-band Hubbard model by Nomura and Yamada.*^^ Ac- 
cording to their third-order perturbation theory, described in 
some detail in § 3, vertex correction terms beyond the second- 
order spin fluctuation terms play the essential roles for the 
pairing; the effects of Coulomb repulsion beyond the boson 
exchange provide the mechanism of unconventional super- 
conductivity. In the heavily doped case as in Sr2Ru04 with 
four electrons with degenerate three bands, a p-wave state is 
more stable than the d-wave states. Odd parity associated with 
the /9-wave symmetry requires the zero-gap at the symme- 
try points at the Brillouin-zone (BZ) boundary; naturally the 
main superconducting gap takes a minimum at the Fermi sur- 
face M-point, because it is closest to the BZ symmetry point. 
The anti-phase anisotropy on the passive bands is also a natu- 
ral consequence if the antiferromagnetic spin fluctuations due 
to the nesting of a and /? bands are not favorable to the pairing. 
Thus, the essential superconducting symmetry is given by the 
chiral p-wave d(k) = zAoik^ + iky), while the low-lying quasi- 
particle excitations are attributable to the nearly-zero gap on 
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Fig. 13. (left) Correspondence between the quasi-particle excitations un- 
der magnetic fields precisely aligned within the basal plane and the mag- 
nitudes of the fourfold oscillation amplitude of the specific heat, (right) 
Dependence on the azimuthal angle (/> of the specific heat at T = 0. 12 K 
for the same reduce field sfi'ength of H/Hc2(T). The disappearance of the 
oscillation at lower fields is attiibutable to the competition between the di- 
rections of the superconducting gap anisotropy in the active (y) and passive 
(a and /?) bands. From Deguchi et al. (2004).*" 



the passive bands forced by the nesting. We will review other 
theoretical developments in § 3. 

2.5 d-vector state 

Theoretically, several contributions to determine the di- 
rection of the (i- vector have been considered. The determining 
mechanisms of the (i-vector orientation can be divided into 
two categories.'*'* One is of the "atomic" origin: orbital wave 
functions of Ru-4c/ and 0-2p affect the anisotropy of the spins 
in both the normal and superconducting states. Such atomic 
spin-orbit interaction (SOI) energy is estimated to be about 
50 meV.'"^ (More recently, the SOI energy is calculated to be 
90 meV at the T point by nonrelativistic band calculations,'"'* 
and to be at least 100 meV based on the comparison between 
de Haas-van Alphen experiments and relativistic band cal- 
culations with local density approximation (LDA).'°^) Such 
SOI energy needs to be considered for electron Cooper pairs, 
while irrelevant to atomic superfluidity of ^^He. In the nor- 
mal state, the anisotropy of the susceptibility is described as 
Xc{T) = xci + l-02;(f„i(r), with the temperature-independent 
part Xi) amounting to about 4% of the total susceptibility.''' 
The anisotropy of the temperature-dependent part, interpreted 
as the anisotropy of the spin susceptibility, is only about 2%. 
Thus in the normal state, there is no strongly preferred orien- 
tation of the spins. In the superconducting state, the lifting of 
the degeneracies of the five possible c/-vector states due to SOI 
decides the preferred orientation of the (i- vector From the mi- 
croscopic theories of the pairing reviewed in § 3, the magni- 
tude of the splitting is expected to be four orders of magnitude 
smaller than SOI itself. 



The other mechanism to align the li-vector is due to dipo- 
lar interaction between the magnetic moments of the spins of 
a Cooper pair in the chiral state, which favors the c/-vector par- 
allel to the L-vector, namely along the c direction. This mech- 
anism dominates the determination of the cZ-vector orientation 
with respect to the orbital moment in chiral superfluid phase 
^He-A. The manipulation of the c/-vector orientation has been 
demonstrated by Ahonen et al}^ For ^He confined within 
narrow parallel plates, the orbital motion is restricted and it is 
expected that the L-vector aligns perpendicular to the plates. 
The fif-vector points in the same direction due to the dipolar 
interaction of about 5 mT. By the application of the in-plane 
field, the c/-vector is indeed reoriented. 

The first mechanism, atomic SOI associated with the Ru- 
4d orbitals, is incorporated in the microscopic pairing mech- 
anism by Yanase and Ogata, '^^^ yielding the c/-vector pointing 
in the c direction with the anisotropy energy corresponding to 
90 mT. Later, Yoshioka and Miyake introduce the importance 
of the 0-2p orbitals in the Hubbard model (the d-p model) 
and reproduced the || c result, but at the same time show 
that d \\ ab becomes more stable if the on-site Coulomb re- 
pulsion at the oxygen site, Upp, is comparable to Udd- A more 
detailed discussion of the atomic SOI will be given in § 3.4. 
The second mechanism, the spin dipolar interaction within 
a Cooper pair, favors d \\ c with the energy corresponding 
to 22 mT. Unlike superfluid ^He, the charges of the Cooper- 
pair electrons lead to the orbital magnetic moment, associated 
with the orbital angular momentum L, which interacts with 
the spins. Such SOI within a Cooper pair yields the spins ori- 
enting parallel to L, thus d \\ ab hy the energy correspond- 
ing to 64 mT. '"^ Within the five c/- vector states allowed in the 
tetragonal crystal symmetry, a ^f-vectorin the basal plane may 
be constructed as a mixing of four of the states. Whether such 
a state is allowed as the ground state is an issue of current de- 
bates. These proposed microscopic mechanisms of orienting 
the (i-vector are summarized in Table V. 

With these competing mechanisms of comparable mag- 
nitudes, it is not easy to predict the c/-vector orientation the- 
oretically. We note that if the c/-vector lies in the basal plane, 
as illustrated in Fig. 1(c), the pinning within the basal plane 
must be extremely weak because of the tetragonal symmetry. 
It would then be possible that the li-vector orients perpendicu- 
lar to the external field of any directions; the spin susceptibil- 
ity is unchanged across T^, consistent with the NMR Knight- 
shift observations for both H \\ ab and H \\ c. 

Collective excitations of Cooper pairs manifest the avail- 
able internal spin and orbital degrees of freedom and give in- 
valuable information of the pairing state. A curious anomaly 
of the spin-lattice relaxation rate of '^O under zero magnetic 
field (NQR) is present in the in-plane oxygen site. Since the 
principal axis for the oxygen quadrupole moment lies in the 
basal plane, l/Ti consists of the spin fluctuations of both the 
in-plane and out-of -plane components. By comparing with 
the ""Ru NQR relaxation rate, Mukuda et al. deduced the 
anisotropy in l/Ti of '^O NQR, which exhibits a broad peak 
below Tc.'"^ Miyake recently"^^ gives explanation of this be- 
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Table V. Various microscopic meclianisms of orienting the rf-vector and the corresponding stability energies. The stability energies are expressed in units 
of magnetic field (mT), divided by the Bohr magneton /jg. *'For Upp > ^tjp. See § 3.4 for details concerning this condition for the oxygen-site Coulomb 
interaction Upp and the hopping parameter tjp . 



Interaction Preferred orientation Stability energy 



Atomic spin-orbit interaction (Ru site) 


d\\L{d II c) 


~ 90 mT 


(Ref. 153) 


Atomic spin-orbit interaction (O site) 


d ±L{d\\ ab) 


~ 50 mT 


(Ref. 160) 


Spin-dipolar interaction within a Cooper pair 


d\\L 


~ 22 mT 


(Ref. 108) 


Spin-orbit interaction within a Cooper pair 


d ±L 


~ 64 mT 


(Ref. 107) 


Zeeman interaction {H \\ c) 


d ±L 






Zeeman interaction {H \\ ab) 


d\\L 







2,0 - • 'J'T.T). 




Fig. 14. Spin-lattice relaxation rate of nuclear quadrupole resonance 
(NQR) of oxygen nuclei. Energy dissipations for both in-plane (open cir- 
cles) and out-of-plane (closed circles) components in these zero-field data 
have been deduced. The anomalous enhancement of the out-of-plane ab- 
sorption is attiibuted to the collective-mode excitation associated with 
the d-vectoT dynamics."" The solid curves are theoretical analyses by 
Miyake.."" 



havior as due to the internal Josephson oscillation of the 
Cooper pairs. As shown in Fig. 14, the oscillation of the d- 
vector in the basal plane, corresponding to the fluctuation of 
Cooper-pair spins along the c axis, leads to l/Ti with charac- 
teristic two broad peaks below T^. Although not available at 
present, '^O NQR signals from the apical oxygen site would 
give confirming results for this interpretation. 

A certain number of theoretical studies on possible col- 
lective modes have been done for the spin-triplet pairing 
state of Sr2Ru04."°""^ Among possible collective modes, 
the most important one is a spin-wave mode which is related 
to fluctuations of the t/-vector direction from its stable direc- 
tion. If the anisotropy in spin space is neglected, then this 
mode corresponds to the Goldstone mode accompanied by the 
symmetry breaking with respect to rotations in spin space and 
its excitation energy equals zero at the long-wavelength limit 
^ = (because the tZ-vector rotates freely without any en- 
ergy cost). In reality, this mode should have an energy gap at 



q - (maybe not zero but small gap for Sr2Ru04), reflecting 
the anisotropy of the c/-vector According to the theoretical 
studies,''"""^ this mode may be observed experimentally in 
the in-plane dynamical spin susceptibility, for the chiral pair- 
ing state d{k) oc z{kjc + iky). Although such a collective spin 
excitation has not been observed experimentally, it might be 
observed in the future, using low-energy neutron scattering 
or other spin-resonance measurements. Discovery of such a 
spin collective mode could be strong evidence for the triplet 
pairing in Sr2Ru04. 

3. Pairing by Coulomb repulsion 

3.1 Microscopic pairing mechanisms based on Coulomb re- 
pulsion 

From theoretical points of view, many studies on the 
triplet pairing mechanism in Sr2Ru04 have been performed 
so far.'^"' '^^ To our knowledge, all of them essentially re- 
quire unconventional (i.e., non-phonon) mechanisms for ex- 
plaining the triplet pairing. It is believed that strong elec- 
tron correlations are much likely to prohibit conventional 
phonon-mediated pairing in Sr2Ru04. In fact, Sr2Ru04 is 
a typical strongly correlated electron system, as suggested 
from many experimental and theoretical studies. According 
to band-structure calculations, '^^ Ku-Ade electrons, which 
are relatively localized, occupy the states near the Fermi level 
dominantly, and therefore the electronic properties will be af- 
fected significantly by electron correlations. Consistent with 
this, sizable mass enhancement was observed by quantum os- 
cillation measurements.'"^'^ In addition, it was confirmed that 
a related compound Ca2Ru04 is a Mott insulator.''*'- '^^ 

Here we overview various proposals on the unconven- 
tional pairing mechanism for Sr2Ru04. Most of proposed 
mechanisms can be classified, depending on the micro- 
scopic origin of the pairing interaction assumed: (i) ferro- 
magnetic spin fluctuations or paramagnons as in the su- 
perfluid 3jjg i 17-121, 137 (-jj-) anisotropic incommensurate an- 
tiferromagnetic spin-fluctuations, '^^"'^^ (iii) incommensurate 
charge or orbital fluctuations,'^* (iv) Hund's-rule coupling 
among Ru-4^fe electrons,'^'- '^^ (v) inter-site interaction (in- 
cluding inter-site Coulomb repulsion), '^^' '^° and (vi) on-site 
Coulomb repulsion between KvL-Ade electrons (particularly 
among 'Ru-Adxy electrons). '■^'"'■^^ At a first glance, (i) may 
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seem the most plausible, since Sr2Ru04 has Fermi liquid pa- 
rameters similar to those of superfluid ^He and the three- 
dimensional SrRuOa is a ferromagnet. However, neutron 
scattering experiments revealed that the dominant magnetic 
correlation in Sr2Ru04 is incommensurate antiferromagnetic 
rather than ferromagnetic,''*^ reflecting the nesting character 
between the a and /3 Fermi surfaces."^"''*'* In (ii), this in- 
commensurate fluctuation is assumed to be much anisotropic: 
Xc{Q)IXah{Q) > 4 - 7 (g is the incommensurate nesting vec- 
tor). '^'^ Qualitatively, such anisotropy is indeed consistent 
with an NMR measurement^^ and theoretical calculation,''*^ 
but seems too large (Xc{Q)IXab{Q) ~ 3 in the NMR exper- 
iment). In (iii), incommensurate charge fluctuations due to 
the nesting between the a and yS bands mediate the triplet 
pairing. In this mechanism, inter-orbital Coulomb interac- 
tion should be larger than the intra-orbital Coulomb repul- 
sion somehow. '^^ In the scenarios (ii) and (iii), the pairing 
interaction will become the strongest on the a and fi bands 
rather than on the main band y, and therefore the a and [5 
bands will be active ones in pairing transition rather than the 
y band, contradicting several experimental results which sug- 
gest the dominance of the y band. In (iv), the Hund's-rule 
coupling is believed to stabilize the parallel spins of triplet 
Cooper pairs. '^^' '^** However, the Hund's coupling between 
the local Ru-4c/ orbitals can stabilize parallel spin states only 
on each local Ru site, but does not necessarily stabilize them 
over the superconducting coherence length. It is still unclear 
whether the Hund's-rule coupling stabilizes spin-triplet pair- 
ing rather than spin-singlet pairing in Sr2Ru04 or not. In (v), 
as easily shown, the in-plane inter-site attractive interaction 
between nearest-neighbor sites induces sin and sin ^^-wave 
pairing. '^^ As for the gap structure in (v), more inter-site cou- 
plings between farther sites become necessary to explain the 
higher-order harmonics of superconducting gap function. If 
the inter- layer interactions between two sites on adjacent dif- 
ferent RuOa layers are included, then horizontal line node ap- 
pears in the superconducting gap.'^''' '■"'■ '''^^ '''^ However, such 
a horizontal node seems inconsistent with a result of field- 
oriented specific heat measurement." 

Recently, Raghu et al. presented a two-band calculation 
for the one-dimensional yz and xz bands, in order to explain 
the absence of edge currents.^" For pairing interaction, they 
used the expression of the second-order perturbation in the 
on-site intra- and inter-orbital repulsion. This theory seems 
substantially the same as in (iii). In fact, the obtained momen- 
tum dependence of the gap function is very similar to that 
obtained in the random phase approximation (RR\) calcu- 
lation of Ref. 126, i.e., the momentum dependence is well 
approximated by the forms A^zik) ~ Aq sin cos ky and 
Av,(^) a! Ao cos A:^ sin A:,, in the both theories. The both the- 
ories suggest that the triplet pairing becomes more favorable 
against the singlet pairings, as the inter-orbital repulsion be- 
comes stronger. 

At present, we consider that (vi) is the most plausible sce- 
nario, in that it adopts only simple assumptions whose origins 
are microscopically clear or which are consistent with most of 



experiments. In the next section we will describe this mecha- 
nism in detail. 

3.2 Single-band Hubbard model 

To describe the above mentioned mechanism (vi), we 
use weak-coupling approximation, but require nothing except 
a realistic tight-binding electronic structure and the on-site 
Coulomb repulsion at each Ru site. First, we consider only 
the y band for simplicity, since the y band takes the major 
part of the density of states (about 57%) and the heaviest ef- 
fective mass (14.6 times the bare-electron mass) at the Fermi 
level. ''"^ We describe the y band by a single-band t-t' Hubbard 
model on a two-dimensional square lattice:'^' 

n.n. n.n.n. 

// = f ^ al^ajo- + t' ^ «L"./(r + fl,'|fl'|flafliT' (4) 

{',j),cr {i,j).cr i 

where a,,^ (a. ) is the annihilation (creation) operator of an 
electron with spin cr at Ru site ;. Here t and t' are the hopping 
integrals between nearest and next-nearest neighbor sites re- 
spectively, U is the on-site Coulomb integral. We take t - - \, 
f - -0.375, which can reproduce the electronic structure 
of the y band approximately. Roughly speaking, since the y 
band is constructed mainly from the K\x-Ad^y orbitals, the mi- 
croscopic origin of U is the Coulomb repulsion between the 
Kw-Adxy electrons, t is estimated to be about 0.25 eV, if the 
bare bandwidth equals 2 eV. "n.n." ("n.n.n.") means that the 
summation with respect to / and j is restricted to nearest (next- 
nearest) neighbor sites. The superconducting transition point 
is determined by solving the Eliashberg equation:'^' 

^<r,a,{k)^-^ ^ V^,^,,^,^,{k,k')\G(k'tA^,^,{k'\ (5) 

A''a"3a"4 

where k - [k, ioj), and Aa-a-' (k) is the superconducting or- 
der parameter (in other words, the wave function of Cooper 
pairs). G{k) is the particle propagator for quasi-particles. 
V,r,(T2.(T4(T,ik,k') is the pairing interaction, in other words, the 
scattering amplitude between pairing quasi-particles, where 
the momenta and spins of the quasi-particles in the initial (fi- 
nal) state are k'crj and -^'0-4 {kcri and -ko-2). Transition 
temperature is the highest temperature below which a non- 
trivial solution Aua-'ik) + appears, and the pairing symme- 
try is determined by the momentum and spin dependences of 
the solution Aa-cr'{k). 

In the weak-coupling regime, we can evaluate the pairing 
interaction by perturbation expansion in U : 

V(k,k') = UV'^^\k,k') + U-V-\k,k') + U^V'^^\k,k') + ■■■ . 

(6) 

Here we have omitted the spin indices for simplicity. The first- 
order term remains non-zero only in the cases of singlet pair- 
ing. This first-order term is independent of momentum and 
frequency, and is repulsive {UV^^\k,k') = t/ > 0) in the s- 
wave channel. This is the reason why the i-wave pairing is 
generally prohibited in strongly correlated electron systems. 
Higher-order terms can possess non-trivial momentum de- 
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pendences, as a result from many-body correlation. The pos- 
sibility of anisotropic (i.e., non-i-wave) pairing depends on 
whether this non-trivial momentum dependence is attractive 
for the corresponding anisotropic pairing channel or not. Thus 
we turn our attention to the higher-order terms. 

The second-order contribution is expressed using the 
bare susceptibiUty xoiq) as U^V^^\k,k') = ±U^Xoik - k'), 
where the sign is plus (minus) for singlet (triplet) pairing. 
In many cases, the bare susceptibility xoiq) captures qualita- 
tively the correct momentum dependence of more precise sus- 
ceptibility xi<l) calculated within advanced approximations, 
e.g., RPA,'"*^ fluctuation-exchange approximation (FLEX),'"^^ 
etc., although the momentum and frequency dependences of 
the accurate xi4) become more enhanced and steeper near 
magnetic transition points than those of the bare xoi^)- I" 
pairing theories based on spin-fluctuation mediation, it is as- 
sumed that the pairing interaction can be expressed by the 
form V{k, k') oc x{k - k'), where xiq) is calculated by those 
approximate methods. '^"^ Therefore, as far as we discuss 
qualitative features including the most favorable pairing sym- 
metry, the second-order theory can often predict qualitatively 
the same results as other methods based on spin-fluctuation 
mediation. These theories present a clear physical picture 
that a pair of quasi-particles are formed into a Cooper pair 
by exchanging elementary excitation ('spin fluctuation' in 
the present case), where V{k,k') is a function only of the 
momentum-energy transfer ^ = k-k', i.e., the momentum and 
energy of the elementary excitation. Pair-scattering processes 
due to exchange of elementary excitations are generally rep- 
resented by such a type of diagrams as in Fig. 15(a). 

Next we proceed to the third-order theory. The third- 
order term U^V^^\k,k') can bring about qualitatively new 
momentum dependence, as mentioned below. The third- 
order term consists of RPA-like terms and so-called 
"vertex-corrected" terms: U^V^^\k,k') = U^V^'^^^\k,k') + 
fy3y(v.c.)|-^^ A:') (Analytic expressions of the third-order terms 
are presented in Ref. 131). The RPA-like terms, which are in- 
cluded also within RPA, present qualitatively the same mo- 
mentum dependence as the second-order term, and are ex- 
pressed by diagrams of the type (a) of Fig. 15. On the other 
hand, the vertex-corrected terms do not depend only on ^ - A:', 
but on both k and k'. Typical examples of third-order vertex- 
corrected terms are shown in Fig. 15(b). To be regarded as 
mediated by some kinds of elementary excitations of the sys- 
tem, V(k, k') must be a function only of the momentum and 
energy [q = k-k') of elementary excitations excited in the in- 
termediate states. Therefore the vertex-corrected terms cannot 
be regarded as due to exchange processes of any elementary 
excitations. 

Within the third-order perturbation theory, Nomura and 
Yamada showed that the most probable pairing symmetry 
on the y band (Ru-4iijy band) is the spin-triplet /9-wave.'^' 
In Fig. 16, transition temperatures calculated for singlet and 
triplet pairing states are displayed as a function of U. For the 
total electron filling of the y band « = 1.33, the highest tran- 
sition temperature is given by a triplet pairing state, while 




Fig. 15. Diagrammatic representations of typical pair-scattering processes. 
The filled four-leg vertices and oriented solid lines denote the on-site re- 
pulsion U and particle propagator, respectively, (a) General diagrammatic 
representation of the pair scattering processes which can be regarded as 
due to exchange of elementary excitations. The rectangle represents arbi- 
fi'ary diagrammatic structure connecting upper and lower parts. Diagrams 
of this type provide contributions depending only on momentum-energy 
transfer q = k-k' for the pairing interaction, (b) Typical third-order vertex- 
coiTected terms. Note that these processes are not included in the processes 
of the above type (a). The diagrams in (b) play an essential role for triplet 
p-wave pairing. 



the singlet Jj.2_j,2-wave state becomes favorable toward half 
filling. Thus, the fact that the electron filling of the y band 
is away from half filling is essential for realization of triplet 
pairing in Sr2Ru04. The momentum dependence of the or- 
der parameter A{k) suggests the p-wave symmetry, as shown 
in Fig. 17 (Only the A:v-type solution At(k) is displayed. The 
A;,,-type solution Ay{k) is obtained by n/2 rotation around the 
z axis. A^{k) and Ay(A:) transform as kj^ and ky under the D41, 
symmetry operations, respectively). Here note that the mo- 
mentum dependence of A^{k) is not similar to the lowest-order 
harmonics sin k^, but seems to contain higher-order harmon- 
ics. Although the dominant harmonics is sin kj^ cos ky, the zero 
contour of A^{k) crosses the Fermi circle only at two Fermi 
surface positions kx - 0. This suggests that line nodes oc- 
cur at kx = on the realistic cylindrical Fermi surface. How- 
ever, it should be noted that these line nodes do not cause 
zero-gap in the energy spectrum, if we also consider Ay{k) of 
the ky type as well, because the real energy gap is given by 
|A(A:)| = [Ax{k)~ + Ay(k)^]^/\^o within weak-coupling theo- 
ries. 

In Fig. 18, V{k,k') calculated within the third-order per- 
turbation theory in U is plotted as a function of k for fixed 
k' . We pay attention to the behavior of V{k, k') near the Fermi 
surface. We should note that V{k, k') takes low values around 
k-k' and the maximum value around k - -k'. The low 
values of V{k,k') around ^: = is attributed to the weak 
ferromagnetic feature of xo(l) (remember that the second- 
order contribution is expressed by U^V^^\k,k') = -U~XQ{k - 
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Fig. 16. (Color online) Transition temperature as a function of the on- 
site repulsion U , calculated by using a single-band Hubbard model for the 
yband.'^' « denotes the total electron filling (« = 1.33 is the realistic value 
for the y band of Sr2Ru04). Half filling is h = 1. 
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Fig. 18. (Color online) fc-dependence of the triplet-pairing interaction 
V{k,k'). k' is fixed at the point indicated by the arrow. The green circle 
represents the y Fermi surface. U = 4 and T = 0.002. 




Fig. 17. (Color online) Momentum dependence of the triplet superconduct- 
ing order parameter A(k) (only the ij-type). The green circle and gray solid 
lines on the base plane represent the y Fermi surface and the zeros of Aijc), 
respectively. 



k')). This ferromagnetic feature is due to the fact that the 
Fermi level is close to the van Hove singularity on the y 
band.'"*^ On the other hand, the high values around k - -k' 
are brought about by the third-order vertex-corrected term 
f;3y(v.c.)(^^ As inferred from the minus sign of the right- 
hand side of the Eliashberg equation (5), this characteristic 
momentum dependence of V(k, k') favors sign change of A{k) 
between k and -k, and accordingly causes p-wave pairing. It 
is interesting that the same vertex-corrected terms as shown 
in Fig. 15(b) cause p-wave pairing also in two-dimensional 
isotropic repulsive fermion systems in the weak-coupling 
regime.'^'' '^^ At present, we do not know how to draw such 
a clear and intuitive physical picture on the momentum de- 
pendence derived from vertex-corrected terms as that on the 
momentum dependence due to elementary-excitation medi- 
ation. The third-order perturbation calculation has been ap- 
plied to other unconventional superconductors than Sr2Ru04, 



high-rc cuprates, organics, several heavy-fermion systems, 
and suggested a plausible pairing symmetry for each of them, 
so far.'^^* Effects of higher-order contributions are discussed 
in the Appendix of Ref. 153. 

Honerkamp and Salmhofer found the p-wave supercon- 
ducting phase by applying the one-loop renormalization- 
group (RG) theory to a similar Hubbard model. '^"^ '^^ Al- 
though the RG analysis does not include one-particle self- 
energy corrections, it has an advantage that the magnetic, pair- 
ing, and other instabilities can be treated on an equal footing. 
According to their result, the momentum dependence of the 
effective interaction seems similar to the one obtained within 
the above third-order perturbation theory, rather than to the 
one induced by the ferromagnetic spin fluctuation, in spite 
that the p-wave superconducting phase is adjacent to a ferro- 
magnetic phase. In fact, the RG calculation seems to include 
the same vertex corrections as the third-order perturbation cal- 
culation. 

5.5 Extension to three-band model 

Nomura and Yamada extended the above single-band cal- 
culation to a three-band calculation.'^^' They used a three- 
band Hubbard model (without SOI) to describe the realis- 
tic three-band electronic structure of Sr2Ru04. In the three- 
band calculation, three components of order parameter, ls.„{k) 
(a = a,/3, y), are required for describing pairing on each band. 
They determined Aa{k) by solving the Eliashberg equation 
within the third-order perturbation in the on-site Coulomb 
integrals. According to their result, the superconducting or- 
der parameter takes the maximum value on the y band. Thus 
the 'orbital dependent superconductivity' initially suggested 
in Ref. 93 is verified by the microscopic theory. The proba- 
ble pairing symmetry is reasonably the p-wave again, and the 
momentum dependence of Ay{k) is almost the same as that 
by the above single-band calculation for the y band. Accord- 
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ing to their three-band calculation, the picture that inter-band 
proximity '"^^ induces pairing on the passive a and /3 bands is 
indeed valid. In fact, they found that the pairing amplitude 
Ay{k) on the y band remains finite and A,j,^(fc) on the passive 
bands become zero, when the inter-orbital Coulomb interac- 
tion is artificially set to zero. 

Nomura and Yamada used the three-band solution Ai,{k) 
{a - a,p,y) for the gap function. "^^^ '^^ If any of the five d- 
vector states displayed in Table IV is realized, the gap mag- 
nitude on each band is calculated by |Aa(^)| = [A^ ^(A:)" -i- 
Aa,y(^)^]'^^U=o, where AaAk) and A.a.y{k) are the kx- and ky- 
type solutions on band a, respectively. The gap function pre- 
dicted by the third-order perturbation theory is presented in 
Fig. 19. The gap magnitude takes the maximum value on the 
y band and the minimum on the fi band, reflecting the domi- 
nance of the y band. The calculated result predicts that the su- 
perconducting gap should be strongly anisotropic in the basal 
plane: gap minimum on the y Fermi surface should exist near 
the zone boundary, and nearly-zero gap minima on the a and 
jS Fermi surfaces should exist at the diagonal points. The gap 
minimum on the y band occurs intrinsically near the zone 
boundary points 0) and (0, since A(^) = holds at 
the points (kx, ky) - (0, 0), (+7r, 0), (0, ±7t) and (+7r, ±7t) in the 
Brillouin zone, as straightforwardly verified from odd-parity 
symmetry and 27r-periodicity of A(A:).'^'' '^^ We should note 
that the nearly-zero gap minima on the passive a and /3 sur- 
faces are not a direct result from the symmetry property, in 
other words, not due to the sign change of the gap function, in 
contrast to the nodes of dj:2_y2-wave pairing. The reason why 
such a strongly anisotropic gap structure appears on the pas- 
sive a and /3 bands is that the antiferromagnetic fluctuation 
due to the nesting between the a and /3 Fermi surfaces weak- 
ens the triplet pairing at these points. Overall, the anisotropy 
of the main gap structure is quite different from that of the 
dxi^yi symmetry 

Nomura and Yamada calculated the specific heat using 
the above gap function. They assumed that the momen- 
tum dependence of superconducting gap near is preserved 
down to low temperatures, and simply determine the abso- 
lute magnitude of the gap below by the standard BCS gap 
equation. The calculated result is shown in Fig. 20, and is 
compared with experimental results. The jump of C/T at 
is dominated by the contribution from the y band. Note 
that, owing to the strongly anisotropic gap structure, the T- 
linear behavior at low temperatures is reproduced even if the 
chiral symmetry is assumed. At low temperatures, thermal ex- 
citations on the a and /3 bands, dominantly contribute to the 
specific heat, due to the nearly-zero gap minima on these pas- 
sive bands. The same gap function can explain comprehen- 
sively the temperature dependences of NMR relaxation rate 
(l/riT),'^^ ultrasound attenuation rate, and thermal conduc- 
tivity.'^''^ 

3.4 Role of spin-orbit interaction and d-vector orientation 

Electron-electron Coulomb interaction is symmetric un- 
der rotations in spin space. Therefore, pairing interaction 



(a) 




(b) 




Fig. 19. (Color online) Predicted superconducting gap structure. • 
(a) Gap magnitude on each of the thi'ee Fermi surfaces is depicted as a 
function of the azimuthal angle </> around the c axis (the a axis direction 
corresponds to </> = 0). The unit of energy on the vertical axis is about 
0.057 eV. (b) The Fermi surfaces are represented by the solid black circles. 
The dependence of gap magnitude on in-plane direction is expressed by 
the distance from the Fermi circle along the direction. 



V(k,k') induced only by the Coulomb repulsion does not 
bring about any anisotropy in spin space, and never lifts the 
degeneracy of the li-vector states in principle. In § 2.5, three 
major contributions to the anisotropic energy for the t/-vector 
have been discussed, as summarized in Table V. Among them, 
we give here a detailed discussion on a leading mechanism 
of lifting the degeneracy, the "atomic" SOI among the Ru- 
4d orbitals. This interaction hybridizes different spin and or- 
bital states with each other and causes anisotropy in spin 
space. Although the SOI is believed to be much 
weak compared with the bandwidth and the Coulomb interac- 
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Fig. 20. (Color online) Specific heat C divided by temperature T is dis- 
played as a function of temperature (normalized by the normal-state 
value). The solid curves and the diamonds (o) represent the theoretical 
and experimental data, respectively. Each contribution of the three bands 
is also depicted separately. 



tion for the Ru-4^/ states (of the order of 100 meV), it plays the 
essential role for determining the most stable cZ-vector state. 
The degenerated J- vector states generally split into five states 
due to the SOI, as shown in Table IV. 

Here we review several microscopic theories including 
SOI. In the early stage, Ng and Sigrist'^^ showed that, un- 
der reasonable assumptions, the SOI among the Rn-A-de or- 
bitals stabilizes the chiral pairing state, where they used sin 
and sin ky basis functions for the gap function. They used a 
realistic three-band model including nearest-neighbor ferro- 
magnetic exchange. According to their results, if the pairing 
interaction is dominant on the y band, then the chiral state 
d{k) = |Ao(sinfcv ± isinA:,,) is reasonably the most stable 
against the others. 

To investigate the possibility of the chiral pairing state 
in a more microscopic level, Yanase and Ogata included the 
SOI to perturbation calculation for a three-band Hubbard 
model. They treated the SOI within second-order pertur- 
bation. For pairing interaction, they performed perturbation 
expansion in the intra-orbital and inter-orbital Coulomb inte- 
grals to third and second order, respectively. The origin of the 
attractive interaction will be the same as that by Nomura and 
Yamada, i.e., the third-order vertex-corrected terms with re- 
spect to the intra-orbital Coulomb repulsion U. The calculated 
eigenvalues of the linearized Eliashberg equation are shown in 
Fig. 21, where note that the c/- vector state presenting the larger 
eigenvalue is the more favorable state. As seen in Fig. 21, the 
chiral pairing state d{k) oc z{kx + iky) is expected to be real- 
ized at Tc- This is an ESP state with the pair spins in the basal 
plane. 

One might naively consider that the magnetic field neces- 
sary to change the direction of the ^f-vector is estimated from 
the SOI energy A via a simple relation ijbB ~ A. Such esti- 
mation is not correct because the anisotropy energy should be 
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Fig. 2 1 . (Color online) Eigenvalues /le as a function of temperature for the 
five rf-vector states:'"^ (1) d(k) oc jtky ± yky, (2) d(k) oc xky ± yk^, (3) 
d(k) oc z(k_y ± iky). The inset shows the same in the wide temperature 
region. 



evaluated from the difference in condensation energy among 
states with different ^f-vector orientations. This energy should 
be of the order of k^ATc, where the difference in transition 
temperature AT^c turns out to be a function of A/W rather than 
A itself, as mentioned below. 

Yanase and Ogata showed that realistic value of spin- 
orbit coupling (about 50 meV) lifts the degeneracy only 
slightly. This is because the first-order contribution with re- 
spect to the spin-orbit coupling vanishes and the symmetry- 
breaking energy is of the order of UJA^/W^ (U: intra-orbital 
Coulomb repulsion, /: Hund's rule coupling, and W: band- 
width), if the hybridization between different orbitals is com- 
pletely neglected. According to their estimation, the splitting 
of Tc between different (i-vector states is estimated as AT^c ~ 
Q.QA-Tc- A very crude estimation using /ibB ~ k^T^- (ATc/Tc) 
leads to 90 mT. Such slight lifting of degeneracy means that 
the li-vector is not so strongly pinned along the most sta- 
ble direction and can be rotated easily even by a small mag- 
netic field. In fact, the J- vector direction must be rotated even 
by a small magnetic field, to be qualitatively consistent with 
the experimental fact that the NMR Knight shift remains un- 
changed at under a small magnetic field (20 mT) along the 
c axis.^-' 

More recently, Yoshioka and Miyake investigated the 
anisotropy of the c/- vector on the basis of the d-p model, in- 
cluding not only the on-site Coulomb repulsion Udd and SOI 
at each Ru site but the on-site Coulomb repulsion Upp at each 
O site.'^" They calculated the pairing interaction within the 
perturbation theory with respect to the on-site Coulomb re- 
pulsions Udd and Upp and the SOI A. They found that, if Upp 
is assumed to be larger than Udd, the c/-vector parallel to the 
ab plane can be more stable than the one parallel to the c 
axis, as shown in the phase diagram of Fig. 22. According 
to their estimation, the splitting of Tc between the d \\ c and 
d ± c states due to the SOI is AT^ ~ 0.022rc for A = 0.2 
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Fig. 22. (Color online) Phase diagram of anisotropy of the c/-vector. The 
stable direction of the rf-vector is parallel (perpendicular) to the c axis in 
the blue (red) region. The green solid line denotes Ujj = Upp points. 



and Udd - Upp — 4.0 (in units of tdp, where tdp is the hop- 
ping between nearest-neighbor Ru and O sites). In addition, 
they showed that the momentum dependence of the order pa- 
rameter changes from sin ky- cos fcy-Hke to simple sin A:^-Uke, 
around Upp ~ 3tdp, as Upp is increased. However, we should 
note that, according to Table IV, any of c/-vector states with 
d ± cdo not break TRS in the absence of magnetic fields, and 
therefore seem incompatible with [iSR results.''^ 

At present, the mechanism based on the Coulomb repul- 
sion at Ru site (the main band is y), incorporated with the 
atomic SOI at Ru site, reasonably supports the realization 
of the chiral pairing state. To be consistent with the NMR 
Knight shift measurement, it is necessary to assume that the 
(i-vector along the c axis should be rotated to the ab plane 
even under a small magnetic field of as small as 20 mT and the 
Coulomb interaction at oxygen site should not be so large as 
the J-vector is laid along the ab plane. Of course, the micro- 
scopic determination of the most stable J- vector state remains 
still a subtle issue, because the stabilization energy of the di- 
rection of the 6^-vector is quite small and moreover possibly 
competing, as mentioned above. 

We would like to comment on one of the promising fu- 
ture directions of the microscopic mechanism theory. In all 
the theories discussed here, simple two-dimensional tight- 
binding electronic structures are adopted. To proceed to a 
more precise evaluation of the mechanism and c/- vector states, 
one should use more precise electronic structures based on 
first-principle calculations. Extension to three-dimensional 
electronic structures may also prove crucial to resolve some 
of the important issues, such as suppression of Hc2 and super- 
conducting multiphase, discussed in § 4. 

Recent LDA calculations of the electronic structure and 
ARPES measurements'"'*' suggest that the three cylindrical 
Fermi surfaces cross near the (110) plane in the absence of the 
SOI, while actually the crossing is removed by the SOI. Ac- 
cording to these results, the orbital character on the Fermi sur- 
face varies with fc;, with a tendency of crossing between two- 
dimensional and one-dimensional bands near k~ = 0. Thus 
we can no longer consider that the y Fermi surface is con- 



Special Topics 

structed almost entirely from the xy electrons. Influences of 
such Fermi surface crossing have not been studied so far, and 
remain an interesting issue to be investigated in detail. 

4. Challenges to the spin-triplet scenario 

In this section, we discuss some of the unresolved issues 
for which explanations in terms of the present spin-triplet p- 
wave scenario have not been settled. We discuss available the- 
oretical models. 

4.1 Suppression of the upper critical field 

The Hc2 of Sr2Ru04 is strongly anisotropic, reflecting the 
quasi-two-dimensional feature of the Fermi surface: 1 .5 T for 
H \\ ab and 0.075 T for H \\ c. The in-plane Hc2 anisotropy ap- 
pears below around 1 K and increases up to 0.04 T. One of the 
unresolved issues on Sr2Ru04 is the temperature dependence 
of //c2 for the field in the ab plane. If the orbital depairing ef- 
fect dominates in determining //c2, Hc2 is expected to increase 
linearly with the initial slope at Tc on cooling and saturates 
to -QJTcdHc2/dT\T=T^ at low temperatures. This behavior 
is described by the so-called Werthamer-Helfand-Hohenberg 
(WHH) theory'^^- '^^ and its extension to p-wave supercon- 
ductors. ' '^^ Although Hc2 of Sr2Ru04 for H \\ c exhibits 
conventional behavior (//c2(0) = -0.78rcd//c2/dr|r=7 ), the 
in-plane Hc2 is strongly suppressed at low temperatures 
(Hc2(0) = -0A2T,dHa/dT\T=Ty' 

Recently, the Pauli paramagnetic effect was suggested as 
a possible origin to limit the in-plane //c2 of Sr2Ru04.'^^' '^^ 
Machida and Ichioka calculated Hc2iff) as well as the field 
dependence of the specific heat and magnetization using 
the Eilenberger equation based on a spherical Fermi sur- 
face with the assumption of the paramagnetic depairing ef- 
fect. Although they suggested that their results reproduced 
some experimental results, the Pauli effect occurs only for ei- 
ther the singlet pairing state or triplet pairing states with the 
c/-vector locked in a particular direction in the basal plane, 
which is inconsistent with the results of Knight shift exper- 
iments suggesting the triplet pairing with the tZ-vector state 
pointing perpendicular to the field in any field direction. The 
NMR result strongly indicates that Pauli effect is not oper- 
ative in Sr2Ru04. Qualitatively similar //c2 limiting is also 
known for another spin-triplet superconductor UPtj for H \\ c 
(Fig. 23),'^^''™ whose origin has not been resolved either. 
In the field direction in which a strong Hc2 suppression has 
been observed, H \\ ab for Sr2Ru04 and H \\ c for UPt3, 
the (i-vector is believed to be perpendicular to the field and 
thus the Pauli paramagnetic depairing is not anticipated. Thus, 
the origins of these Hc2 limitings observed in two plausible 
spin-triplet superconductors remain challenging issues to be 
resolved. 

In order to clarify the detailed features of the //c2 lim- 
iting in Sr2Ru04, Kittaka et al. performed ac susceptibility 
measurements with a precise control of the magnetic field di- 
rection using a vector magnet system. They investigated 
the dependence of //c2 on temperature and the angle 6 be- 
tween the magnetic-field direction and the ab plane, as shown 
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Fig. 23. (Color online) Field-temperature phase diagrams of Sr2Ru04 
compared with those of UPtj . Suppression of Hc2 of similar magnitude 
is evident for H \\ ah in Sr2Ru04 and for H \\ cm UPt3. From Deguchi et 
al. (2002, 2004)'"'"' and Dijk etal. (1993). 



in Figs. 24(a) and 24(b). It was revealed that the //c2 limit- 
ing is remarkable only for \9\ < 5° and remains strong even 
at relatively high temperatures close to Tc. In the previous 
study, a clear kink in Hdiff) around 0-2° was reported 
from the specific heat measurements at 0.1 K,'^' but such a 
clear feature was not detected at any temperature either by the 
thermal conductivity'^' or recent ac susceptibility'^"' '^'* (the 
inset of Fig. 24(b)). It was also revealed that the angular de- 
pendence of determined from the ac susceptibility was 
reproduced with an effective-mass model for an anisotropic 
three-dimensional superconductor, though the anisotropic ra- 
tio needs to vary with temperature. '^^^ It is essential to develop 
a new theory which explains the observed //c2 behavior as 
well as other key experimental results, such as invariance of 
the Knight shift. 

4.2 Superconducting multiphase 

A multi-component superconducting order parameter 
leads to an additional phase transition in magnetic fields. 
For example, another spin-triplet superconductor UPt3 is well 
established to have additional superconducting transition in 
magnetic fields, resulting from its multi-component order pa- 
rameter. Also in Sr2Ru04, an additional transition was ob- 
served just below //c2 at low temperatures below 0.8 K 
(Fig. 25) by the specific heat,'^' dc magnetization,'^^ and ac 
susceptibility.'^'*' '^^ Intriguingly, this additional transition as 
well as the in-plane //c2 anisotropy is suppressed by a slight 
misalignment of the field from the ab plane; they are observed 
only for 0< 1°.'^^ 

Although the mechanism of this superconducting multi- 
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Fig. 24. (Color online) (a) Anisotropy of the upper critical fields //c2 of 
Sr2Ru04. h' is H^i normalized by the initial slope in the vicinity of T^. 
Strong suppression of Hd occurs for the polar angle only within a few 
degrees from the basal plane, (b) Field-angle 9 dependence of H^2 nor- 
mahzed by H^iiS = 90°). The //c2 suppression at 6 = 0° is evident even 
at relatively high temperatures. From Kittaka et al. (2009)."^ The inset 
compares the Hc2iO) determined from various experiments. 
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Fig. 25. Field-temperature phase diagram of Sr2Ru04 detaihng the regions 
of superconducting double transitions for field accurately aligned within 
the basal plane. The phase diagram is constructed from the double-peak 
structure of the temperature dependence of the specific heat. From Deguchi 
et al. (2002).'" 
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phase has been actively discussed, it is still controversial. For 
a superconductor with d{k) oz z(kx + iky), if the magnetic field 
is applied parallel to the x axis, the degeneracy between zk^ 
and zky is lifted and the transition from z(fcj ± ie^y) to zk^ oc- 
curs and the non-chiral single-component zkj^ is expected to 
become stable near at all temperatures.'^^ This additional 
transition can be easily suppressed by a sUght misalignment 
of the field of only one degree. However, in this theoretical 
model, large in-plane anisotropy for such additional transition 
is expected, '^^ which is incompatible with the experimen- 
tal results. Another possibility is the crossover from zkx to 
the spin-polarized (z - iy)^.! state. In this theory, the tran- 
sition from z{kx + iky) to zk^ is additionally expected in the 
region far below Hc2- Although a tiny magnetization anomaly 
far below //(;2 wasreportedby Tenyaef a/.,'^^ there is no other 
report which supports the results of the dc magnetization ex- 
periment. At present, this scenario is not well established, ei- 
ther. Machida and Ichioka proposed that the second transition 
is attributed to the FFLO in Pauli limited superconductors.'^^ 
Clearly, there remains a contradiction between the PauU ef- 
fect scenario and the results of Knight shift experiments, as 
written in § 4. 1 . 

5. Novel superconducting phenomena 

5.7 Superconductivity in eutectic systems 

The intrinsic Tc of pure Sr2Ru04 is well established to 
be 1.5 K.^- However, a large diamagnetic signal and the rapid 
decrease of resistivity with an onset of nearly 3 K are ob- 
served in some single crystals of Sr2Ru04. Soon after the dis- 
covery of such "3-K phase", it was noticed that the enhanced- 
Tc superconductivity originates in the Sr2Ru04-Ru eutectic 
solidification system. '^"^ 

The Sr2Ru04-Ru eutectic system is obtained during the 
single crystal growth of Sr2Ru04 with a floating zone tech- 
nique using Ru self-flux.'^' Because volatile Ru02 evaporates 
from the surface of melt during the growth, excess Ru needed 
as a flux tends to solidify in the core region of the crystal rod 
and coexists with Sr2Ru04, forming a periodic pattern charac- 
teristic of eutectic solidification. Figures 26(a) and 26(b) rep- 
resent polarized optical microscope images of polished sur- 
faces of a eutectic crystal parallel and perpendicular to Ru 
lamellae, respectively.'^" A number of Ru lamellae (bright 
area) with typical dimensions of 30 x 10 x 1 ^m-* are em- 
bedded in a Sr2Ru04 single crystal (dark area). The distance 
between the nearest lamellae is on typically about 10 (im. The 
orientation of Ru lamellae with respect to the crystal axes of 
Sr2Ru04 varies among different crystals and even within the 
same crystalline rod. Surprisingly, although the values of the 
unit-cell parameters are different between Sr2Ru04 and Ru, it 
was revealed from transmission electron microscopy (TEM) 
images that the interface between Ru and Sr2Ru04 is atomi- 
cally sharp, as shown in Figs. 27(a) and 27(b). '^^ By contrast, 
dislocations caused by the lattice mismatch are distributed 
on a larger scale, as displayed in TEM pictures as dark lines 
(Figs. 27(c) and 27(d)). 

Various experimental results suggest that enhanced-Tc 



Fig. 26. (Color online) Optical microscope pictures of the surface of a 
Sr2Ru04-Ru eutectic crystal. Thin Ru metal lamella are embedded in the 
matrices of Sr2Ru04 single crystal. From Maeno et al. (1998).'*" 




Fig. 27. Transmission electron microscopy images of Sr2Ru04-Ru eutec- 
tic crystals. From Ying et al. (2009).'*" Note that the dislocation lines in 
(d) are the results of superposition over many layers. 
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superconductivity is non-bulk and non-i-wave, occurring in 
the Sr2Ru04 region in the eutectic system. First, a broad 
transition with a large diamagnetic shielding'**" and a tiny 
anomaly of the specific heat around 3 K"*^ indicate that the 3- 
K superconductivity is non-bulk. Second, a zero-bias conduc- 
tance peak (ZBCP) observed above the bulk of Sr2Ru04 
in the Sr2Ru04-Ru junctions suggests that non-i-wave su- 
perconductivity occurs in the 3-K phase. "^"^'^^ Third, the 
anisotropy of Hd clearly depends on the orientation of the 
crystal axes of Sr2Ru04 rather than that of Ru lamellae,"*'^' '^^ 
indicating that Sr2Ru04 is a stage of enhanced-Tc supercon- 
ductivity. In the 3-K phase, the formation of Josephson net- 
work consisting of inter-lamellar supercurrents is suggested 
by the zero-voltage current observed at temperatures well 
above the bulk of Sr2Ru04"*'* and a large diamagnetic 
shielding with a characteristic response to ac and dc magnetic 
fields.''*" 

One of the peculiar features of the 3-K phase is tem- 
perature dependence of H^a-^'^^ Using a definition of //c2 as 
the resistivity inflection point associated with the supercon- 
ducting transition to the 3-K phase, //c2 appears to behave 
as (1 - T/Tcf^^ for H || c and (1 - T/T^f-'^^ for H || ab 
close to Tc.'^^ In addition, HdiT) for H \\ c was revealed 
to have an upward curvature around 2 K. Such unusual be- 
havior of Hc2 is explained using a phenomenological theory 
within the Ginzburg-Landau formalism for an interfacial su- 
perconductivity by Sigrist and Monien (S-M).'"' They as- 
sumed one interface and analyzed with the vector order pa- 
rameter d(k) = zirjxkx + Tjyky), where (z/j., 77,,) = ( 1 , i) is plau- 
sible in Sr2Ru04. In the S-M model, it was suggested that the 
presence of interface parallel to the y axis makes one of the 
order parameters ky more stable than the other ^i. By intro- 
ducing a finite width of the interfacial superconductivity in 
the S-M model, Matsumoto et al. succeeded in reproducing 
the unusual temperature dependence of H^^i more quantita- 
tively.'''^ However, the mechanism of the strong in-plane H^^i 
suppression at low temperatures and the origin of the hys- 
teresis behavior observed just below H^2 for H \\ ab remain 
unresolved."*^ 

By analyzing the interplay between the chiral p-wave or- 
der parameter of Sr2Ru04 and the i-wave order parameter of 
a cylindrical Ru inclusion based on the Ginzburg-Landau for- 
mulation, Kaneyasu and Sigrist pointed out that, at a tempera- 
ture somewhat lower than of Ru (~ 0.49 K), a state emerges 
with a spontaneous flux distribution around a Ru inclusion.'''^ 
Moreover, Kaneyasu et al. proposed that a discontinuous 
transition from time-reversal-symmetry-conserving phase to 
time-reversal-symmetry-breaking phase occurs on cooling at 
a temperature between 1.5 and 3 K.'"'' Their model explains 
anomalous asymmetric features with respect to the current di- 
rection in the current-voltage characteristics observed below 
2.3 K'**** as well as the observation of ZBCP below 2.4 K."*^ 

Recently, hints for clarifying the mechanism of the en- 
hancement of were provided by Kittaka et al. from uniax- 
ial pressure experiments. '''^ They found that uniaxial pressure 
along the c axis of only 0.2 GPa induces a broad supercon- 
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Fig. 28. (Color online) Diamagnetic shielding signals of a pure Sr2Ru04 
crystal under uniaxial pressure. Evidence for the enhancement of supercon- 
ducting transition temperature up to 3 K is clearly visible. From Kittaka et 
al. (2010)."^ 



ducting transition with the onset of about 3.2 K even in a 
pure, Ru-inclusion-free Sr2Ru04 sample with a sharp transi- 
tion at Tc - 1.35 K at ambient pressure. The results of the 
ac susceptibility measurements are shown in Fig. 28. It was 
also found that both out-of-plane and in-plane uniaxial pres- 
sures largely enhance the magnetic shielding fraction associ- 
ated with 3-K superconductivity in eutectic samples. '"^''"^ In 
particular, in-plane uniaxial pressure is much more effective 
to the enhancement. By contrast, hydrostatic pressure sup- 
presses the bulk Tc of Sr2Ru04'''^' '"" as well as the shield- 
ing fraction associated with 3-K superconductivity.-"*^ These 
facts strongly indicate that an anisotropic crystal distortion of 
Sr2Ru04 enhances its T^. Nevertheless, the microscopic ori- 
gin of this enhancement by anisotropic distortion has not been 
clarified. 

Let us describe another Sr2Ru04-based eutectic system, 
Sr2Ru04-Sr3Ru207. Fittipaldi et al. found that Sr2Ru04- 
Sr3Ru207 eutectic crystals can be grown with a feed rod con- 
taining 45% excess Ru with respect to Sr2Ru04 under a gas 
mixture of 10% oxygen and 90% argon with the total pres- 
sure of 10 atm.^"' X-ray diffraction analyses of Sr2Ru04- 
Sr3Ru207 eutectic crystals indicate that the directions not 
only of the c axis but also of the in-plane axes of Sr2Ru04 
and Sr3Ru207 are common in the eutectic system.'"' It was 
also revealed from TEM images that interfaces parallel to the 
growth direction (longitudinal direction of a crystalline rod) 
are sharp and defect-free whereas interfaces perpendicular to 
the growth direction appear wavy and decorated with Ru pre- 
cipitates of micron-size.^"^ In addition, mono-layers of RUO2 
planes with a thickness on the order of nm (very thin Sr2Ru04 
layer) were revealed to be embedded as stacking faults in the 
apparent Sr3Ru207 regions of the eutectic system. Neverthe- 
less, the Sr3Ru207 domains of the eutectic system have a sig- 
nificantly lower level of impurities, such as Sr4Ru30io (n - 4) 
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and SrRuOs (« = oo), compared to SraRuaO? single-phase 
crystals. ^'^■^ 

In ac susceptibility measurements with an apparent 
Sr3Ru207 domain cut from the Sr2Ru04-Sr3Ru207 eutectic 
system, multiple superconducting transitions at 1.3 and 1 K 
were observed. ^"^'^"^ These transitions were revealed to be 
non-bulk from specific heat measurements and are consid- 
ered to be attributed to the superconducting transitions of thin 
Sr2Ru04 layers embedded in the Sr3Ru207 domain as stack- 
ing faults. In fact, the anisotropy of //(;2 of these transitions is 
the same as that of the bulk Sr2Ru04 and H^i for H \\ c of 
the transition at 1.3 K (1 K) is compatible with (sUghtly lower 
than) that of the bulk Sr2Ru04.^"'^'^°'' 

At the early stage, the possibility of the long-range prox- 
imity effect in the eutectic system arising from the bulk 
Sr2Ru04 region into Sr3Ru207 was discussed because finite 
critical current was observed in a crystal containing a large 
amount of Sr3Ru207.^°^ This finite critical current is proba- 
bly not due to the long-range proximity effect, but due to the 
presence of very thin Sr2Ru04-layer inclusions. In fact, 
of the multiple superconducting transitions is not very sensi- 
tive to the electronic state of Sr3Ru207; although Sr3Ru207 
exhibits ferromagnetism under the uniaxial pressure along 
the c axis above 0.3 GPa,^"**'^"^ the onset temperature of the 
superconducting transition at 1 K is almost invariant or in- 
creases slightly under the uniaxial pressure up to 0.4 GPa.^'° 
This uniaxial pressure effect is consistent with the scenario 
that the multiple transitions originate from the embedded thin 
Sr2Ru04 layers. 

Recently, Yanase focused on a model of a thin Sr2Ru04 
layer intercalated in the Sr3Ru207 domain as a multi- 
component spin-triplet superconductor with the random spin- 
orbit coupling.^" He proposed that the li-vector state in such 
a Sr2Ru04 layer is different from that in the bulk Sr2Ru04 
and suggested that the cZ-vector state is locked to be d{k) oc 
xky — yky- with time-reversal symmetry. 

5.2 SriRuOix as a topological superconductor 

Recently, quantum phenomena for which underlying 
topology of wave functions plays fundamental roles are ac- 
tively investigated with a universal point of view. Topology is 
a concept and method enabling the classification of shape by 
continuous deformation. For quantum phenomena, the shape 
subject to continuous deformation is based on the symme- 
tries and phases of the wave functions. A well-known ex- 
ample of topological quantum phenomena is vortex quanti- 
zation in a superconductor, for which phase winding associ- 
ated with real-space topology determines the state. In addi- 
tion, novel quantum phenomena may emerge in materials for 
which bulk wave functions themselves have nontrivial topol- 
ogy, characterized by winding numbers defined in momentum 
space; thus these materials are called the "topological mate- 
rials". Topological materials include some insulators, super- 
conductors and superfluids and are classified in terms of dis- 
creet symmetries of the Hamiltonian as well as the special 
dimensions.^'^ It has been recognized that states specific to 
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non-trivial topology and immune to disorder emerge at the 
edge, i.e., its surface or interface with other materials, accord- 
ing to "bulk-edge" correspondence.^''' 

A chiral p-wave superconductor with the wave function 
kx + \ky is classified among such topological materials and is 
called a topological superconductor, along with other classes 
of topological superconductors. For unitary-spin, chiral fcv + 
\ky state with broken time-reversal symmetry, the topologi- 
cal invariance is specified by the Chern number.^"* The chiral 
edge state discussed in § 2.3 is a consequence of bulk-edge 
correspondence reflecting the bulk topological properties of 
a chiral p-wave superconductor. The existence of such chiral 
edge sate has been demonstrated by a recent experiment of 
quasi-particle tunneling spectroscopy.^'** In superconductor- 
insulator-normal metal junctions with Sr2Ru04 and Au as 
shown in Fig. 29, broad humps filling up the gap are ob- 
served in the conductance spectra. Compared with the narrow 
ZBCP associated with the flat-band edge state of YBa2Cu307 
at its (110) interface, the broad conductance peak is a direct 
manifestation of the chiral edge state characteristic of chi- 
ral p-wave superconductivity as shown in Fig. 30. In addi- 
tion to a number of evidence indicating chiral p-wave sym- 
metry, this confirmation of the chiral edge states reassures 
that Sr2Ru04 is classified as a topological superconductor. 
Very recently, Sasaki et al. found that a bulk superconduc- 
tor Cuj:Bi2Se3,^'^'^'^ derived from a doped topological insu- 
lator Bi2Se3, exhibits a similar ZBCP in a "soft" point-contact 
junction with gold wire attached on the cleaved surface of a 
Cu vBi2Se3 crystal by silver paste. -'^ After examination of the 
tunneling spectra in terms of possible superconducting states 
based on the relatively simple band structure expected for 
Cuj^Bi2Se3, they concluded that the observed ZBCP gives cru- 
cial evidence for topological superconductivity with odd par- 
ity in Cu vBi2Se3. 

Such ZBCP emerging at a junction between an unconven- 
tional superconductor and a normal metal is a consequence of 
the density of states due to Andreev bound state. ^"^ It has been 
recognized that such Andreev bound state and the penetration 
of the pairing amplitude into the normal metal by the proxim- 
ity effect are intimately related. For an unconventional spin- 
singlet superconductor, the formation of the Andreev bound 
state at the junction surface of the unconventional supercon- 
ductor and the proximity effect are competitive, but it turns 
out that for a spin-triplet superconductor, they may be con- 
structive, making the proximity effect unusually strong and 
long-ranged.^'^ As one consequence of such enhancement, it 
is predicted that due to the proximity penetration of the An- 
dreev bound state into the normal metal, in a T-shaped junc- 
tion consisting of a spin-triplet superconductor and a normal 
metal, the density of state of the normal metal increases below 
of the superconductor. 

Moreover, it has been recognized that such Andreev 
bound state can be described in terms of odd-frequency pair- 
ing states. In a diffusive normal metal, states other than 
isotropic i-wave states become difficult to survive. At a junc- 
tion with a spin-triplet superconductor, therefore, spin-triplet 
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Fig. 29. (Color online) (a) Scanning ion microscopy image of a Sr2Ru04- 
Au junction fabricated on a Si02 substrate and a schematic of the junction 
structure, (b) Comparison between experimental conductance specti'a at 
0.55 K (solid lines) and theoretical spectra (dotted lines). From Kashiwaya 
etal. (201 1).^"* 



i-wave, namely an odd frequency pairing state, is realized. 
Intimate relations among surface bound states, proximity 
effects, odd-frequency pairing states, and topological edge 
states are discussed in a recent review by Tanaka, Sato and 
Nagaosa.^'-' 

Concerning the superconducting junctions, it has been 
recognized that small-size crystal is important to extract the 
intrinsic unconventional behavior of Sr2Ru04. For a focused 
ion beam (FIB) fabricated small crystal of Sr2Ru04/Ru eutec- 
tic, an unusual hysteresis in the current- voltage characteristics 
has been observed. It was interpreted as due to the motion of 
the chiral domain walls which alters the critical current. 

As a novel phenomenon ascribable to real-spatial topol- 
ogy of chiral superconductivity, the critical current of a junc- 
tion between Sr2Ru04 and the 5-wave superconductor Pb 
has been investigated. Pb/Ru/Sr2Ru04 proximity-junctions 
(Fig. 31) utilizing Ru inclusions in the Sr2Ru04/Ru eutectic 
crystals as a means of achieving a good electrical contact be- 
tween the two superconductors were fabricated. In such junc- 
tions, unusual temperature dependence of the critical current 
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Fig. 30. (Color online) Angle-resolved conductance spectra for (a) d^y- 
wave superconductors and (b) chiral p-wave superconductors. The ver- 
tical axis is the quasiparticle energy normalized by Ao. The brightness 
represents the magnitude of the junction conductance. Bright regions cor- 
respond to the conductance peak originating from Andreev bound states 
formed at the edge. Conductance spectra for (c) c(vy-wave superconductors 
and (d) chiral p-wave superconductors. From Kashiwaya et al. (201 1).^'^ 



is obtained: it increases on cooling from about 3 Kbut exhibits 
an unusual sharp drop at the onset of bulk superconductivity 
of Sr2Ru04 at 1 .5 K; at lower temperatures the critical current 
increases again. This unusual behavior is ascribable to topo- 
logical competition of the phase winding numbers around a 
Ru inclusion of proximity-induce i-wave superconductivity 
in Ru and of superconductivity in Sr2Ru04 as illustrated in 
Fig. 32; thus the junction may be called a topological super- 
conducting junction. The phase winding number of Sr2Ru04 
changes from that of non-chiral p-wave superconductivity of 
the 3-K phase to the chiral p-wave superconductivity of the 
bulk Sr2Ru04. We note that a similar phenomenon was re- 
ported before in Pb/Sr2Ru04/Pb junctions with the critical 
current measured between the two Pb electrodes.^-'' In the re- 
cent investigation, it is demonstrated that this behavior is es- 
sentially obtained in a Pb/Ru/Sr2Ru04 junction configuration 
with one Pb electrode. 

Another remarkable phenomenon associated with real- 
space topology of superconductivity in Sr2Ru04 is emergence 
of a state related to half-fluxoid state originating from the spin 
and orbital degrees of freedom. We will describe it in the next 
section. 

5.3 Half-quantum vortices 

One of the truly remarkable phenomena expected for 
spin-triplet states is emergence of the half-quantum vortices 
(HQV) associated with the magnetic flux just half of the flux 
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Fig. 31. (Color online) Variation of the critical current 4 with temperature 
of a Pb/Ru/Sr2Ru04 junction. The inset is a schematic of the junction, 
/c sharply drops just below 7'(;,sr2Ru04 but increases again below a certain 
temperature designated as T' . From T. Nakamura et al. (201 1).~^^ 
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Fig. 32. (Color online) Schematic images of the evolution of the order pa- 
rameter at the Sr2Ru04/Ru interface. In the upper panels, the gradations 
in the small circles depict the momentum-directional dependence of the 
superconducting phase (p at each spatial position; the arrows denote the 
momentum direction for which = 0. The angle is defined as nor- 
mal to the interface (see (a)). The lower panels represent the supercon- 
ducting phases (p{8) at the Sr2Ru04/Ru interface under no external cur- 
rent. The solid lines represent the phase of .v-wave superconductivity in 
Ru and the broken lines that of p-wave superconductivity in Sr2Ru04. (a) 
Tc,Sr2Ru04 <T <Ti,: The irj ± iekj^ state with the winding number W = is 
realized (A' phase), matching with the i-wave order parameter induced in 
Ru. (b) T < 7'c,Sr2Ru04 : Replacement by kx±\ky, the bulk state of Sr2Ru04, 
with W = ±1 (B phase), (c) T <«. rc,sr2Ru04: Increasing interfacial energy 
enlarges the phase deformation in the .v-wave, strengthening the Josephson 
coupling. From T. Nakamura et al. (201 1).^^^ 



quantum Oq = hjle. It was proposed in the context of su- 
perfluidity of 3]^g 227,228 i^^j j^^g never been confirmed 
experimentally. In general, flux quantization occurs as a 
consequence of the single-valuedness of the order parame- 
ter around a singularity; either a hole in a superconductor or 
the vortex core created in the superconductor. Ordinarily, this 
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single-valuedness is satisfied if the orbital-phase (f variation 
around a singularity, 6ip = j> V^-ds, is equal to 2nn with an in- 
teger n. However, for a spin-triplet superconductor of the ESP 
type, the single-valuedness of the order parameter may also be 
achieved by a combination of the variations of the spin part as 
well as the orbital part. To demonstrate this fact, we consider 
an ESP spin-triplet state |*I')esp with the spin quantization axis 
along the x axis, which can be expressed as 



l*I'>ESP = V2^e'^ (0, cosQ', sin a) , 



(7) 



where the (i-vector direction a corresponds to the phase of 
the spin part as represented in eq. (A. 2) in the Appendix. The 
single-valuedness of this state is satisfied not only if 6(f - 2nn 
and Sa - Q but also if 6ip - {2n - l)n and Sa - {2m - 1)jt, 
where m is an integer The latter is an HQV, because only 
the orbital phase if couples with the vector potential and thus 
the magnetic flux. The simplest example with n - m - I is 
illustrated in Fig. 33(b), in which the (i-vector is represented 
by an arrow. As we shall see in the Appendix, |*I')esp can also 
be expressed as 



l*l')ESP=^(-e''MTT).v + e'«MU).v), 



(8) 



where 6-^ - (p - a and 0| - ip + a. Thus, an HQV can be inter- 
preted as a vortex for which winding numbers of the compo- 
nents |TT).i and |J,J,),i are different. Specifically, an HQV with 
the minimal vorticity, namely with |2« - 1| = \2m - 1| = 1, 
should be coreless in the sense that one of the two compo- 
nents |tT)x and lii) V does not have singularity at the center of 
the HQV, as illustrated in Fig. 33(c). 

The HQV involves the rotation of the d-vectoT and the 
associated spin currents expressed as oc Vo-, in which a 
is the directional angle of the cZ-vector Thus the total energy 
associated with an HQV is 



-EhQV = C[-pcharge(VvP)^ + Pspin(VQ;)^] 



(9) 



with a coefficient C. The charge current associated with an 
HQV is screened in a distance characterized by the London 
penetration depth, whereas there is no such screening effect 
expected for the spin current. To limit the increase of the en- 
ergy due to spin current, it is anticipated that either a meso- 
scopic sample size or a closely-separated pair of HQVs is nec- 
essary. ^^"^ In the former, when the size of a sample is compa- 
rable to the penetration depth, HQVs are similar in energy to 
full-quantum vortices (FQVs). 

In search of an HQV, measurements of magnetic torque 
of micron-size rings of Sr2Ru04 have recently been per- 
formed. As shown in Fig. 34, a micrometer-size crystal with 
a hole drilled by an FIB technique was placed on a silicon 
cantilever of laser-detected ultra-sensitive torque magnetome- 
ter^' The purpose of the central hole is to avoid the complica- 
tions associated with vortex core production. In such a small 
sample, the charge supercurrent is non-zero everywhere and 
it is a fluxoid, given by the sum of flux and kinematic term, 
rather than a flux, that is quantized to an integer multiple of 
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Fig. 33. (Color online) Illustrations of the phases of the superconducting 
wave function around (a) a full-quantum vortex (FQV) and (b,c) a half- 
quantized vortex (HQV). As illustrated in (b), a ;r-rotation of the J-vector 
direction around an HQV and the w-shift of the orbital wave function re- 
cover the single-valuedness of the total wave function. The panel (b) dis- 
plays an example of an HQV in which the phases of the order parameter 
change from (if, a)=(0, njT) to (n, "injT). In this example, d rotates from 
d II z to II -z- An HQV can be also interpreted as a vortex for which 
winding numbers of the components |TT>.v and lii)j are different. An ex- 
ample equivalent to the case in (b) is shown in (c). Here, 6| = - a is 
constant whereas 6; = i/j + a changes from ;r/2 to 5;r/2. Note that the 
actual phase for |tT).v is 0^ + ;r due to the minus sign in the first term in 
eq. (8). 




Fig. 34. (Color online) Silicon cantilever used in an ultra-high-sensitivity 
torque magnetometer with a micron-size crystal ring of Sr2Ru04 glued. 
The Sr2Ru04 ring is also shown in the inset. From Jang et al. (201 1).'' 



flux quantum: 



As - n<l>, 







(10) 



With increasing the magnetic field penetrating the hole, the 
magnetic moment of the ring is found to decrease uniformly 
with a negative slope corresponding to the Meissner screen- 
ing. In addition, a positive jump is observed periodically, cor- 
responding to the full fluxoid quantization. The interval of 
these jumps is consistent with the fluxoid quantization given 
by eq. (10). 
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Fig. 35. (Color online) (A) Magnetization of a micron-size ring of a 
Sr2Ru04 crystal under magnetic fields H-, showing a slope due to Meiss- 
ner current as well as exhibiting jumps associated with fluxoid quantiza- 
tion. With additional fields perpendicular to the ring axis, the fluxoid jumps 
split into two. (B) Magnetization of the ring sample after the Meissner 
shielding contribution has been subtracted. The magnitude of magnetiza- 
tion steps of split jumps is nearly half of that corresponding to a full flux- 
oid jump. The plateau thus is consistent with a formation of a half-fluxoid 
state, closely related to tiie HQV state. From Jang et al. (2010).^' 



Unusual jumps between the FQV jumps are observed 
when additional field is applied perpendicular to the axis of 
a hole as in Fig. 35. The range of longitudinal field for the 
stability of the half-fluxoid regions increases nearly propor- 
tionally to the perpendicular field. Experimentally, the size 
of each magnetization jump as well as the magnitude of the 
longitudinal field is precisely consistent with the half fluxoid 
expectation. Other possible reasons of the fractional magne- 
tization jumps were examined. For example, penetration of 
quantized vortices through the side walls may result in the 
magnetization jumps of varied magnitude. To test this possi- 
bility, the perpendicular field is applied to different in-plane 
direction on the same crystal, but the character of the magne- 
tization jumps remains unchanged: it always occurs as half- 
fluxoid jumps. The measurements were repeated in several 
crystals: the HQV feature was not detected in a crystal with 
the size greater than a few microns, or in NbSe2 crystals with 
a comparable anisotropy as Sr2Ru04. 

The origin of the necessity of the perpendicular fields 
has not been fully understood. A semi-quantitatively feasible 
mechanism is due to "kinematic polarization" proposed by 
Vakaryuk and Leggett.^^' In an HQV state, the orbital wind- 
ing numbers of the spin-up and spin-down species differ by 
unity. The associated difference in the kinetic energy due to 
unbalanced charge current results in the unequal population 
of each spin species. Application of the field parallel to the 
ab plane would then lower the spin-polarization energy of an 
HQV state compared with that of an FQV state. In order to es- 
tablish the HQV state, other direct measurements of the flux 
quantization are much desired. In any case, the observations 
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of half-height flux jumps by Jang et al. certainly give addi- 
tional strong confirmation of the spin-triplet ESP state real- 
ized in Sr2Ru04. 

6. Establishing spin-triplet superconductivity 

We have reviewed the current status of our knowledge 
of unconventional superconductivity in Sr2Ru04. From the 
material perspective, the investigations are greatly helped by 
the availability of large size, high-quality single crystals: they 
are chemically-stable, stoichiometric, and having the quasi- 
particle mean free path more than an order of magnitude 
greater than the superconducting coherence length. The su- 
perconducting eutectic systems containing Sr2Ru04 add op- 
portunities to study novel superconducting properties includ- 
ing exotic proximity effects. From the viewpoint of the elec- 
tronic structure, strongly-correlated multiband Fermi surfaces 
are relatively simple thanks to the structural simplicity, and 
the normal state properties are quantitatively well character- 
ized as a Fermi liquid. This allows development of micro- 
scopic mechanism theories of unconventional superconduc- 
tivity based on the realistic multiband Fermi surfaces. 

From the extensive NMR studies, the realization of the 
spin-triplet pairing seems quite certain. Nevertheless, obtain- 
ing the complete picture of superconductivity of Sr2Ru04 still 
remains as a challenge. As a promising future direction, we 
have emphasized the necessity of including spin-orbit interac- 
tions in conjunction with three dimensional aspect of the elec- 
tronic states. Inclusion of these factors may lead to improved 
interpretations of the mechanism of superconductivity as well 
as crucial understanding of some of the unresolved issues. 
To advance further, it is necessary to revisit the supercon- 
ducting double transitions near H^i as well as the collective- 
mode excitations observed in the NMR relaxation under zero- 
magnetic field. Relatively unexplored areas also include the 
properties of vortex lattice especially with the field parallel to 
the basal plane. 

A new perspective as a "topological superconductor" 
should also be important. Generation of half-quantum vor- 
tices, observation of phenomena associated with chiral edge 
states, fabrication of superconducting junctions exhibiting 
novel behavior such as odd-frequency pairing states are some 
of the promising future directions. Studies using micron-size 
crystals are proven to be essential to uncover some of these 
phenomena. At the same time, in order to firmly estabUsh the 
chiral p-wave state with the two-component order parameter, 
superconducting properties under symmetry-breaking fields, 
for example, in-plane uniaxial pressure, may provide useful 
information. 

The extensive efforts to elucidate the superconductivity 
of Sr2Ru04 in all details demonstrate the possible approaches 
that can be taken to unveil novel behaviors of an unconven- 
tional superconductor In this sense, we hope that this review 
also serves useful for the studies of other unconventional su- 
perconductors. 
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Appendix: Two equivalent interpretations of the half- 
integer quantum vortex 

As briefly mentioned in § 5.3, there are two equivalent in- 
terpretations of the half-integer quantum vortex (HQV): One 
is that (2« -1)71 phase shift in the orbital part of the order pa- 
rameter and additional (2m - \)n phase shift in the spin part 
around an HQV. The other interpretation is a vortex for which 
winding numbers of the two ESP components are different. In 
this Appendix, we explain the equivalence of these two inter- 
pretations. 

First, we write again the definitions of the bases of the 
d-vector: 

i = |5, = 0> = ^(-|TT>z + IU).), 

V2 

y = |5v = 0) = ^(ITT>. + IU>.), 
V2 

z = |5, = 0> = ^(in), + liT>,), 
V2 

where the subscript z in the state vectors indicates that the 
spin quantization axis is the z axis. If one takes x as the spin 
quantization axis, an alternative equivalent set is obtained as 

i = |5, = 0> = ^(in). + liT>.), 

V2 

5) = i5v = o> = ^(-inxv + iuxv), 

V2 

z = = 0> = ^(ITT>.v + IU).v). 
V2 

Here, we start with a spin-triplet equal-spin pairing (ESP) 
state |*I')esp with the quantization axis along x: 

I^)esp = A (-e'^MTT). + e^Mii).) ■ (A-1) 

We note that the minus sign in the first term is just for a con- 
venience that the state with the zero phases 0-\ - Oi - Q cor- 
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responds to the y state. Then, we can convert them to another 
set of phases: a = {0i - 0-\)l2 and ip = {0i + 0^)/2. Using a 
and (fi, |*F)esp can be expressed as 

I'I')esp = ^e'^ (-e-'"|TT>.v + e'"lii).v) • (A-2) 

This equation can be further modified using the relations |tt 
h = -(y + mi V2 and lii), = (3) - mi V2 as 



I*I')esp = [{e'" + e-^y + {e'" - e-)(-iz)] 
= V2zfe''^(0, cos a, sino-). 



(A-3) 



Thus, this relation indicates that ip, the average of the phases 
of the up- and down-spin components, corresponds to the 
phase of the orbital part of the order parameter, and that a, 
the half-difference between 0-f and Oi, corresponds to the di- 
rection of the (i-vector in the plane perpendicular to the quan- 
tization axis. 

When an ordinary full-quantum vortex (FQV) exists, the 
orbital phase (p must shift by +2«7r after a circulation around 
the vortex in order to satisfy the single-valuedness of the 
wavefunction. According to eq. (A-3), the single-valuedness 
is also satisfied if the phases a and ip both additionally shift 
by +71. This is the HQV, because only the shift in the orbital 
phase (fi couples to the vector potential resulting in the mag- 
netic flux of 0o/2 - hIAe. As one can see, the n phase shift 
in a around an HQV corresponds to a rotation of the li-vector 
by 180° around the quantization axis. These facts are illus- 
trated in Fig. 33(b). In an alternative viewpoint, the phase for 
the up-spin component |TT)a is - f - and the phase for 
\ii)x is 0J, - ip> + a shown in eq. (A-2). Thus, the winding 
numbers of 9-\ and 9i are different around an HQV. A remark- 
able example is the case that the shifts in a and ip around an 
HQV are both +n. In such cases, only one of |TT)a and |iJ,)i 
exhibits a phase shift of +27r, whereas the other exhibits zero 
shift. Thus, one of the two components |TT)a and |ii).v has 
zero vorticity around an HQV as schematically explained in 
Fig. 33(c). Such an HQV is "coreless" in a sense that one of 
the two components has no singularity at the center of the 
HQV, if it is realized in a simply connected bulk sample. 
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